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ABSTRACT 


The replenishment at sea (RAS) maneuver is studied in 
detail fcr heading and speed eentrol. Design Of: 
purposefuliy ncnlinear control laws is accomplished fcr the 
Mariner hull using the linearized equations of amecticn in 
three degrees cf freedon. Extensive use of lew order 
modeling and optimal control theory was made. Procedure 
steps are resented in detail to facilitate redesign for 
other ship tyres. The results are verified using DSL 
Simulaticn fer a number of possible RAS scenarios. The 
ccntrol systems are also tested in a sea state to insure 


prcper cperation in the presence of external perturbations. 
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I. INTRODUCTION 


The advent ef the” ‘digital computer as standard 
equipment cn -toard virtually all modern Naval ships has 
opened the field of Digital Computer Contrcl in almost all 


aspects cf ship life. The computer has been a viable asset 


in fire ccrtrol systems for years and has been used 
extensively for aids to ship maneuvering control in the form 
of NIDS (Naval Tactical Data System) readouts. The 
declining ccsts cf generai and special purpose computers has 


made their inception as a manpower replacemert or 


augmentaticn a reality. Their high speed and accuracy can 
make them perform functions with much greater safety than 
previously attainaktle with time proven (and sometimes time 


weary) “"seaman's eye." 


This then is the basis for this thesis; a study of the 
types of maneuvers that can be handled more accurately and 


safely than presentiy being accomplished. 


One area of study is the total Replenishment At Sea 
(BAS) proklem including both course and speed control for 
tke apprcach and alongside phases. This ‘Sztuaticn. has 
always been cne of extreme danger due tc the collision 
Fotentials invclved. However, other dangers are invclved in 
the on deck evclutions when the ships are not kept at a 
fairly ccnstant distance. Sudden violent maneuvers may 
cause the stress on the lines between ships to increase 
enough t¢ cause the lines to part. The reality of this 
danger is readily apparent to anyone who has ever seen a 
Marila line cr steel cable part or a kingpost shackle break 
or a kingfost suddenly tend under these extreme stresses. A 


system which will minimize these potential dangers is well 
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» | ey mero 
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werth investigation. 


Gf course with a digital computer, the algorithm for 
RAS can ke modified or replaced to enable its use as a 
Maneuvering contrcl device for other situations such as 
formation steaming or single ship navigation transit 


Gentrol. 
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II. MODELING 


A. MARINER CYNAMICS 


In the conception of this thesis, realistic models of 
modern destrcyer hull configurations were sought. This 
search fproved fruitless. The hydrodynamic coefficients for 
present cay destroyers are not currently available. 
However, scme naval and civilian research is presently being 


conducted to cktain these coefficients. 


A cemgplete set of these coefficients is necessary for 
any maneuvering control system design. A hull configuration 
which has teen under continual study with well defined and 
verified hydrodynamic coefficients was chosen{!]., This hull 


is commonly referred to as the "Mariner" hull. 


The development of the equations of motion in six and 
three degrees of freedom have been well documented[?]. The 
mcdel used fcr this thesis is the equations of moticn in 
three degrees of freedom linearized with second order and 
higher terms eliminated. These equations are characterized 
by dependency on small perturbations about a specific 
operating point. The maneuvers experienced in the fcllowing 
chapters dc not entirely meet this criterion. The 
inadequacy and shcrtcomings of this model are of little 
censeguence Lecause relevant hydrodynamic ccefficients are 
not available, and the methods presented can be applied to 
any ship tyre. 


The development of the model is readily available to 


the interested reader [3]; only a summary of the equations 
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and their ccrresfonding hydrodynamic ccefficients 


presented kere. 
The equations of motion used are as follows: 


X.-£ U+X U-u ) +X d=0 

( u ah 9? d 

(Y.-m) V+¥ V#(Y —-m)R+Y.R+Y d=0 
V v E G d 


(No-I )E+N R4#N.V+N V+N ~d=0 
EZ e v Vv d 


The direction and sense of the terms in the 
eguations are shown in figure II-1. Letting: 


a =E-Y. 
11 Vv 
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Figure II-1l 


Direction and Sense of Hydrodynamic Terms 
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The equaticns cf motion can Le written as: 


"t "w 
a Atb Atc Ata Bteb Btce B=IF1 
11 1 11 2b 24 ZA 


“a w 
a Atb Atc Ata Beb Btc B=IF2 
TZ 1z Z Ze 22 22 


a Ctb Cte C=IF3 
33 S13) 


33 
er: 
n " 
a Ata: E=£1 
ZA 
iT i] 
a Ata E=1I2 
V2 2 
" 
&. £C=£3 
gs 
where: 
IF1=-Y ed 
IF2= N ed 
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ze3s- y ode 
: u 


and: 


I1=-E i-c A-b B-c BtIF1 
oth 2H OF 


By sclving this system of equations, the following 
relaticnshifs are established: 


A=(a a_ It-a_a_ 12)/DEL 
Zo. 33 sai 2.4 


B=(a.a I2-a a go ADE 


DEL=é (aa =a vat) 
SSaeolte22* 42> 2 


which yields the solution: 


a Tc. on 
veiea +f A dt 


t, a t z t 
J=B=9 + B dt=9§ [ CB (0) +f Bdt jdt 


° € W 
U=C=u f Clade 
04 


The space cocrdinate system is defined as fcllows: 
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Y=UeSINY+Vecosy 


X=U «COSY-VeSIND 


where: 


Table II-4 summarizes the symbols and ncemenclature used 
in the fcregcing abbreviated solution of motion in three 
degrees of freedom. The applicable characteristics of the 
Mariner hull are presented in table [II-2 with the 
corresponding nondimensionalized hydrodynamic coefficients 
and the DSL computer program variable names delineated in 
Sapte: II-3. 


Computer Program #1 is the basic DSL program that was 
developed frcem these equations of motion. This progragz uses 
two ships to illustrate the turning characteristics cf the 
Mariner hull fer various rudder commands. Figure II-2 shows 
the difference between a step mcdel and a ramp model rudder 
in a gecgraphic flot. Figure II-3 shows the correspcnding 
difference in yaw. 
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TABLE II-1 


SYMBOLS AND NOMENCLATURE 


Definition 


derivative 


respect to 


derivative 


respect to 


derivative 


transverse 


derivative 


transverse 


derivative 


te angular 


derivative 


to angular 


derivative 





of longitudinal force with 


longitudinal acceleration 0 


of longitudinal force 


longitudinal velocity 


ef lateral force 


acceleration V 


oF lateral force 


velocity V 


of lateral force 


acceleration R 


GE Lateral force 


velocity R 


cf lateral force 


to rudder angle d 


with 


with 


with 


with 


with 


with 
U 


respect to 


respect to 


respect 


respect 


respect 


derivative of yaw moment with respect to 


transverse 


acceleration V 


derivative of yaw moment with respect 


te transverse velocity V 
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TABLE. Ei-fF (cont) 


SYMBOLS AND NOMENCLATURE 





SymrEcl Definition 
ie derivative of yaw moment with respect to 
£ 


angular acceleration R 


N derivative of yaw moment with respect to 
angular velocity R 


N derivative of yaw moment with respect to 
rudder angle d 


R yaw angle acceleration 
R yaw angle velocity 
oe initial velocity of origin cf body axes 


relative to: fluid 


V transverse acceleration of ship axes 
relative to flui 





Vv transverse velocity of orign of ship axes 
) relative to fluid 


x hydrodynamic longitudinal fcrce 
yy hydrodynamic lateral force 
0 longitudinal acceleration cf ship axes 


relative to fluid 
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TABLE II~1 (cont.) 


SYMBOLS AND NOMENCLATURE 


Symbcl Definition 





a 


OLE arcs, velocity of ship axes relative 
te fluid 


q yaw angle 

At actual time 

T nondimensionalized time 

= longitudinal distance that the ship center 


of gravity is forward of the ships axes 


u longitudinal velocity of ship axes relative 


to fluid (ceperating point) 


TABLE fI-2 


CRARACTERISEICS OF MARINER-TYPE SIUDY SHIP 


Length,.,£t). 52728 
Beam, ft 716...0 
Drakt,, £¢ 29575 
Displacement, tons 16,800. 
Block coefficient, ce 0.6 


a en ae es ial 


ache «2 owas eA gel 
¢ 


Cav , taeios ak 
d =f 





TABLE II-3 


NCNDIMENSICNAL HYDRODYNAMIC COEFFICIENTS 


Computer Program Nondimensional 
Coeficient Variable Name Value 
(X.-B) MXUD -0.0085 
u 
X XU =0).00 12 
u 
¥ XV. ~0.01243 
v 
(Y.-m) MYVD -0.015 
Vv 
(Y -mu ) MYR -0.0051 
rE 1 
(Y..~-8x. ) YRD —0.00027 
r g 
> YDELR +0.0027 
N NV ~0.00351 
v 
N. NVD -0.000197 
v 
(N -mx u ) NR -0.00227 
© 1 
(veri) IZNRD -0.00068 
Nie 
- NDELR -0.00126 
X SS -0.0000462 
n 
¥ —— =0.0000052 
n 
N --- +0.0000026 
n 
an XDELR 0.0 








TABLE LE=37 (Cont .) 


NCNCIMENSIONAL HYDRODYNAMIC COEFFICIENTS 


Values LEased on the following operating point: 
a =1.0 (15 Kts) 
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BF. RUDDER RESPONSE 


The previous section indicates a marked difference in 
behavicr tetween step and ramp rudder models. This prempted 
an investigation into realistic rudder modeling which would 
fulfill tke reguirements of limit stops and maximum tcudder 


rate. 


NSRCC[*] has modeled the rudder of the DD-931 Class 
Destroyer. The basics of this model are presented in the 


block diagram of figure II-4. 


The first limiter models the rudder stcrfs which fcr the 
Mariner are +30 degrees. The second limiter models the 
proportional band of a variable-displacement pump by 
iworting its gaximum percent stroke. The limits for this 
nenlinear element have been found to be +7 degrees. 

The transfer function ae) accepts an input error 


Signal of up to 7 degrees, converts it to a rudder rate, and 


integrates the rudder rate to obtain rudder angle. Letting: 
d = Maximum rudder rate (2.0 degrees/sec) 
m 


al = Maximum error input (7.0 degrees) 
emax 


The system gain can te defined as: 


K =4d /4a 
g & Emax 


= 2.0/7.0 
= 0.2€5714 /sec 
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Tc convert this model to the required nondimensionalized 
form, the fcllcwing manipulation is required: 


K *' = K eL/u 
g g 1 


" 


0.2857140e527.8/ (15¢1.689) 


5.95224 


it} 


where: 
L = shif length 


ie = operating point speed (15 Ktse1.689 ft/sec/Kt) 


Computer Program #2 is the DSL program which sodels 
this system. The curves of figure II-5 and II-6 exhibit the 
responses cf various step rudder commands. These are 
takulated and cress referenced in table II-4. 


These responses show the characteristics of a realistic 
rudder in that the rudder is never allowed tc slam intc the 
stcps. They exhibit the time delay between command and 
response which is a function of the rate of response (2.0 
degrees/sec). A control system design with this scheme is a 
much more difficult problem than one with an idealized 
Cudder (step response) because the entire rudder ccntrol 


system beccmes quite nonlinear. 
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TABLE II-4 


RUDDER COMMAND AND RESPONSE 


Rudder Ipitial 
Figure Curve Command (deg) Condition(deg) 
Ls 1 +30.0 -30.0 
Tr—5 -2 +#25.0 . -25.0 
II-5 3 beso -20.0 
T—5 4 +15.0 —15.0 
II-5 5 +10.0 -10.0 
TI-5 6 + 5.0 - 5.0 
II-6 1 +30.0 0.0 
II-6 2 +2500 0.0 
II-6 3 #20.0 0.0 
II-6 4 +15.0 0.0 
II-6 5 +10.0 0.0 
II-6 6 + 5.0 0.0 
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ce ENGINE RESEONSE 


Figure II-7 pertrays a complex model of a gas turbine 
pecpulsicn fiant{ Sj[*]. This model contains the elements 
required for a complete dynamic study of the system. For 
the purpese cf this thesis, such a complicated model is not 
required if the overall input-output relationship can be 
estaklisted. 


Reference 5 establishes an - output speed (U) 
relaticnship for a step input of desired speed (Ud) and is 
redrawn as figure II-8. The relationship aprears to be that 
of a first cr seccnd order system with a time delay. 


The system equations for a first order approxisation 
With a time delay may be written as: 


SPLIN = K*SPDDES*e 3 

SELERR = (SFDIN-SPDCUT) *G 

SPLOUT = Jf sepeRe dt 
Which yeilds the transfer function: 


[Ss 
SELTCUT (s) G*K*e 


SPLIN(s) s+G 


paiich is bleck diagrammed in figure II-9. 





Figure II-9 
Propulsion Plant Low Order Model Block Diagram 
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From figure II-8, the time delay, system gain, and time 
ccenstant can be estimated as: 
T = 4.88 sec 
0.9877 
0.092 


Computer Program #3 was used to okttain the ster input 
response. The original complex system output and the low 
order approximaticn are compared in figure II-10. As 
indicated in this figure, the two responses are very close. 
Considering the linearized approximations made in the 
equations cf motion, this response is accurate enough for 
system study use and is used as the model fcr speed ccntrol 
tm Chapter III. 


Similar methods may be used to obtain simplified 
low-order mcdels for other. high-order propulsion systems now 
in use (e.g. pressure fired boiler systems, 1200 1b. 
systems, etc.). They may not, however, sigplify to a first 
order aprfrceximaticn suitable for system study. A method of 
computer determination of low-order models of high-crder 
systems is contained in ref. 7 and may be mechanized for 
this purpese. 
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D. EXTERNAL FORCES 


The modeling of ship dynamics cannct be complete 
without the introduction of external forces which perturb 
its respcnses. These forces are caused by many factors and 
scme are mcre relevant than others in the scope of this 
thesis. The two that are considered can cause substantial 
perturbations that must be modeled and eventually acccunted 
for in tke ccatrol system design. 


1 Iwo Ships in Proximity 


—— eee =-—_—_- — 


Whenever two ships operate in close proximity (less 
than 250 feet), suction and pressure forces between hulls 
are present. Studies have been conducted on the Mariner 
mOlif[*] which have produced data for ccnstruction of a 
farily of curves fcr two ships passing on the same heading. 
No data has keen found for the cases of two ships not cn the 
same heading. Other restrictions on the work presented in 
Ber. tt are that the ships are of the same type and of 
Similar hull ratios. 


Interactive effects vary as the square of speed. 
However, this is only true if both ships are at the same 
speeds. The interaction modification factor is based cn the 
normalizeé speed of 15 kts. This factor can thus be written 
as: 


2 
SPDE = CDOT! 


Exact effects on the interacticn forces and moments in the 
Situation where the ship's speeds differ are not available. 
This is inconvenient since during the approach phase, the 


ncrmalized speed of the approach ship (ship B) can te as 
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high as 1.5. If the effect on ship B is as stated above but 
with its cwn speed causing the interaction modification, the 
interactive forces and moments can be 2.25 times greater 


than withcut sfeed considerations. 


Withcut the ability to pin down this relationship, it 
was ignored in the development cf the contrcl laws presented 
in this thesis. Appendix C was written with the expressed 
mtent of illustrating the effect of modifying the 
interactive forces and moments to the extremum case 
mentioned arkcve. It must be realized that this case is not 
ccnsidered likely in that it is felt that the interactive 
forces and mcments modification cn ship B are more apt to be 
caused Ey the speed of ship A. If this is so, since ship A 
is kept at a constant 15 kts., the interactions need no 
mcdificaticn fcr speed consideration in this thesis. 


Reference 1 also gives a method of modifyin¢e the 
interactive forces and moments based on different ship 
lengths. Fcr ease of computation and graphical 
presentation, the two ships were considered of equal 
lengths. Te modify this to ships of dissimilar lengths, the 
resulting hydrcedynamic derivatives must be mcdified as shown 
in ref. 1 (alsc shown in appendix C). 


Since ne closed form expression existed for these 
forces, tke family of curves reproduced as figures II-11 and 
II-12 were guantized in subroutine SLOPES (an adaptation of 
the subroutine of the same name from ref. tt)... (Appendix °C 
contains a curve fitted subroutine that was compiled after 
tke completicn of the research on this thesis. It was not 
used for any design c@ simulation runs except for those 
presented in that same appendix.) An interpolation 
algorithm is used to approximate the intermediate values 
between quantized values and between the curves of the 


family. 
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* 
Note: 


Lateral Distance DY (Ft) 


-400 -200 0 200 400 


To convert to normalized lateral distance DY - 
Divide by the ship length L. 


Figure II-1ll 


Family of Interactive Y force Curves 
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* 
Lateral Distance DY (Ft) 


“400 200 0 200... ~+~-400 


* 
Note: To convert to normalized lateral distance DY - 
Divide by the ship length L. 





Pagure: ii=12 
Family of Interactive N Moment Curves 
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The main purpcse of this subroutine is to compute the 
interactive forces between ships in the replenishment at sea 
Situation anc output the values for perturbation of the 
ccentrol ship only. Even though both ships are affected by 
‘these perturkations, a one ship control system which is 
effective regardless of the other ships’ moticn is 
considered. Ccnseguently, the interactive forces on the 


seccnd ship are ignored. 


Subroutine SLCPES is contained in appendix A. Figure 
Ii-13 is a geographic plot of the two ships passing at 105.6 
feet with their rudders amidships (0 degrees). Figure II-14 
and figure II-15 show the magnitude of the lateral force (Y¥ 
force) and rctational moment (N moment) of the reference 
ship on the ship making its approach (control ship). The 
reference ship is at 15 kts. and the apprcach ship is at 
Zoat akts'. fhe, Toeontro:, ship ‘starts its approach 5 ship 
lengths (2639.0 feet) astern and 0.4 ship lengths (211.12 
feet) laterally displaced. The most graphic portrait cf the 
effects cf these fcrces and moments appears in the yaw 
changes which are presented in figure II-16. From these 
figures it becomes readily apparent that these perturbations 
cause viclert motions of the ship which must be acccunted 
for in any contrcel system development. Throughout the 
development of stch a control system in chapter III, these 
forces and mucments are considered inherent in the model for 
RAS ccntrcl. 
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2. Waves 


The nodeling of sea state in the form cf waves and wave 
interacticns has cccupied the time of many researchers 
ae The exact formulaticn of waves will act be 
accomplished in this thesis. Since the main concern here 
will be to test the control scheme developed in chapter III, 
a much sigzplified wave generator can be used. To introduce 
the required experimental perturbations cn the designed 
control system a periodic wave system with a fundamental 
frequency ard its second harmonic is used. Some small 
Candom wave properties are introduced that ride on these two 
Sinusoids. A simple expression of this combined wave can be 


written as: 


2 
W = WResin(WE)+(ELlewWF /WL) esin(2eWE) +WFeWRVeSin (WE) 


W denctes the Wave 

WF is the Wave Force 

WE is the Wave Encounter radian Frequency 
WL is the Wave Length 

WRV is the Wave Random Variable 

PT is 2..1415926 


With this wave as a basis, a method cf mcdeling this in 
the dynamic environment of the total ship simulation was 
defined. The modeling includes the intrcduction of this 
wave intc the three degree of freedom equaticns of mction. 
To accomplish this a set of defining relationships were 
develcped. First the general wave direction is input to 
establish the direction of the wave encounter on the ship. 
If the ship direction is YAWDP2 and the wave directicn is 


WD, the expected wave direction is defined as: 


ExCL = WD-YAWDP2 
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Next the wave enccunter frequency (radian frequency) can be 
estaktlished with knowledge of the ships normalized true 
speed (CDCT2), wave length (WL) and normalized wave velocity 
(WV). The wave encounter frequency (WEF) is then: 


WEF = ZePle (CDOT2+WVecos (EWD))/WL 


The total wave encounter (WE) is ncthing more than the 
wave enccunter frequency (WEF) times time. This gives the 
wave enccunter radian frequency required in the simple 
expressicn fcr the combined wave previcusly shown. 


This dces not ccmplete the task, since the individual 
wave forces of each degree of freedom must te derived for 
this general wave expression, namely the ccmponents of WF. 
Again a much simplified version of the more complex real 
life wave forces were used. The  6X+ rand Y¥) -Loreées’ ~are 
considered first. These can ke modeled as cosine and sine 
functicns cf the expected wave direction (EWC) such that: 


WFEX 
WEY 


WFecos (EWD) 
WFesin (EWD) 


where WF is the total wave force of the enccuntered wave. 


The rotational N forces are a little more difficult to 
establish. Ey ccnsidering that no rotational forces are 
created Ey a wave directly on the bow or stern, or directly 
cff the bean, and that it iS maximum when the wave is at 45 
degrees off the bow OE stern, a much simplified 
approximaticn is developed. Realizing that this method is 


very crude, the N force can be written as: 


WEN = WFesin (2eEWD) 


To add mcre creditability to the wave mcdel, a random 


4g 





variable is added to the wave force at the waves fundamental 
freguency. A gaussian (normal) distribution was chosen with 
a zero méear and a standard deviation of one-tenth the 
Maximum allowatle force of WF. A zero mean signifies that 
the expected amplitude of the random wave is 0.0, while the 
standard deviaticn signifies that 68% of the random waves 
will have agplitudes less than one-tenth of the maxinun 
allcwable fcerce of WF. Also, 94% will have amplitudes less 
than one-half the maximum allowable force of WF. This small 
added perturbation allows for verification of the model 
Sexegulaticn with a ‘stochastic force, which in turn adds 
creditabkility to the developed ccntrol systems. 


What remains is to define the total wave force (WF). 
It is impertant not to fall into a common simulation pitfall 
which inevitably causes unneeded design changes. A sea 
state does nct increase at an infinite rate. It therefore 
Someincorrect. to start a simulation with initial conditions 
set for calm sea and immediately introduce a high sea state 
perturbaticn. The initial large perturbaticn transient can 
give results that are not only unrealistic, but can cause 
the model and control system to produ3e unstable results. 
This is especially true in this case since the linearized 
(stall fperturkaticns about an operating pcint) equaticns of 
mcticn are used. 


With this in mind, a ramp feed in of the wave force 
with a limiter at the desired maximum wave force (WFMA) was 
used. The slcpe of the ramp was established to impart 
Minimum initial transients, yet increase the wave force to 
an accertable testing level within the time frames cf the 
Simulaticnos cf chapter III. The slope is designed such that 
the maximum wave force is reached in 94.815 seconds actual 


Simulaticn time (4.548 seconds problem time). 


Computer program #4 was used to verify the wave action 
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model. Table II-5 on page 56 indicates the figures preduced 
and changes in input wave length (WL) and wave direction 
(WD) for e€ach run. The input parameters that were constant 
for all runs are tabulated below: 

YAKDEF2 = 0.0 

CEGT2 = 1.5 


aa 
WS = 5.0 
WFMA = 0.1137 


* 
NCTE: WS is the unnormalized wave speed. Conversicn to 
ncrg@alized wave velocity is: 
WV = WS/15.0 


Intrcduction of the wave forces is accomplished by 
multiplying the established wave forces by the rudder 
hydrodynasic coefficients for the individual reference 
directions. This effectively scales the wave forces tc the 
ship model téing used. The wave force result is coded in 
tke ship simulation program as follows: 


IF12 = KA1eD2+YY2+KATeWY 
IF2Z = KE1eDZ+YN2+KBIOWN 


IF32 = kKC1eD2+NC2+KC1 WX 


Detailed results of the wave farce effects are given in 
Chapter III and will not be dealt with here. 
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* 
NOTE: 





WL is given in ship lengths 

WD is given in degrees 

Curve numbers of all runs corresponding to 
wave force components are: 


Force Curve 
WX al 
WY 2 
WN 3 
Table II-5 


Wave Simulation Listing 
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IIIT. REPLENISHMENT AT SEA 


A. HEADING CCNTROL 


fen cControl Cheice 





Many studies involving replenishment at sea (RAS) have 
treated the preblem asa multivariable system[{ 19 j[11}[12] 
{13]. Academically, there is nothing wreng with thus 
approach. Hcwever, aS a practical system it leaves much to 
be desired. The key drawback in the multivariable system is 
the inescapaktle dependency on a command and control link 
between tke replenishment ships. The unreliability of UHF 
communications at these close distances is a much 
experienced phenonenom. It is felt that any knowledgeable 
ccmmanding cfficer would not entrust the safety of his ship 
to such a guestionable link. An alternate method which is 
described kere is a modern extension of the long trusted 
“seaman's eye" ccncept, where the senscrs and ccntrol 
devices umtst be self contained on the Ship making the 
Bpproach (tereafter referred tc as the Eeceiving Ship oL© 
ship B). 


In all present day RAS operations, the ship on which 
the apprcach is heing made (hereafter referred tc as the 
Supplying shir or ship A), must maintain the replerishment 
course and speed. The receiving ship accomplishes the 
Maneuvers tc maintain station relative tc the supplying 
ship. 


The parameters which are presently measured "visually" 


are relative position (in both the X and the Y directions), 
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relative heac (usually in reference to crdered replenishment 
course), and relative motion in the X direction (fcr speed 
matching). These parameters are usually visualized by the 
ccnning cfficer who in turn gives corrective orders to the 
helmsman. The helgsman must then translate these verbal 
orders into rudder and speed ccmmands through the helm and 
lee helm ccnscles. The accuracy of the execution of the 
conning cfficer's crders is extremely dependent cn the 
ability cf the helmsman and throttleman. This system can be 
quite effective, and it can also be quite disastrous. This 
fast reacting and constantly changing environment lends 
itself tc kEreakdcwn in communications and manifests the 
inability of scme individuals toa cope with the required 


critical ran-sachine interfaces. 


To eliminate these problems, present state of the art 
digital ccmputers and sensors are available for immediate 
implementaticn of a completely automatic ship Gentrol 
systen. Such a control system may be installed on 
individual ships and be used for RAS without the reguirement 
of having the matching installation on the cther ship cf the 
heckup (another drawback of the multivariable approach). 


One cf the many pitfalls that may be encountered in 
digital simulation is the reality of the parameters that are 
measurable ir the real world situation as opposed to those 
that are incidently available in the simulation. With this 
fact as a keynote, Subroutine RBMEAS (Range and Bearing 
MEASurement) was developed. This subroutine, as listed in 
appendix A, defines the forward (FWD) and after (AFT) 
relative and true bearings, and ranges from the receiving 
Soap tO Similar points on the supplying ship. Figure III-~1 


delineates tke terms used in the subroutine. The SDFn terms 
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Measurement Techniques 
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are used tc position the bow and stern sensors and 
reflectors cn ship B and A, respectively, in geographic 
cocrdinates as a function of ships head. The FWD distance 
on the X cccrdinate is ADFX and the Y cocrdinate is ADFY. 
Sizilarly, tke AFT distances are ADAX and ADAY. R1 and R2 
are the FWE and AFT ranges measured by a highly accurate 
Canging device installed on ship 8B. This same tfranging 
device, if properly provided with a pinpoint reflector on 
the supply ship(ship A), will give accurate relative 
bearings FWD and AFT., B1 and B2 respectively. The distance 
between senscrs may be varied, but as a rule should te kept 
as far apart as pessible to allow maximum sensitivity. The 
distance used in this thesis is 1.0 (one ship length), and 
the distances were considered the same for beth ships. This 
is net a necessary condition and may be changed to suit the 


situaticn. 


Subreutine RBFEAS assumes highly accurate sensors in 
both rance and Eearing measuring ability. Such sensors are 
presently available in the form cf Radar altimeters{1*] and 
Laser ranging devices. Another possibility for a measuring 
method is a single sensor time sharing to chtain range and 
Heating tc kEcth reflectors from a single device. Such a 


Sirgle senscr scheme is sketched in figure III-2. 


Once tke FWD and AFT parameters are available, they may 
then be used to determine other desired quantities. 
Subroutine HDGRAS (HeaDinG control for RAS) was develofed to 
Output tke desired heading corrected for heading difference 
G@eeeship A and’ B and the prejected correcticn for distance 
€rrer. This subroutine is listed in appendix A. The center 
Tange and tearing are the average of the FWD and AFT range 
and bearing cutput from Subroutine RBMEAS. 


The additional heading due to distance is projected as 


if ship E maintained its present course until it was 
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Figure III-2 
Alternate Method of Measurement 
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perpendicular to the center of ship A. The reasoning behind 
this is illustrated in figure III-3. If the present course 
will cause ship B to arrive on the station desired (DS), no 
heading chance is required. The expression for PSIACC (J 
Additioral heading due to Distance Correction) is: 


PSIADC = RSENS (DDC+DAeSIN (AA1) ) 


where: 
RSENS = Range SENSitivity gain 
DDC = Distance Desired Corrected for side of apprcach 
DA = center Distance Absolute (range) 


AA1 = 360 degrees - relative bearing of center position 


The headince difference of ship A and 8 is desired 
Since, even if the range when alongside is correct, a large 
disparity in heading cannot be tolerated. It is realized 
that some heading difference (crabbing) is necessary to 
Maintain tke distance. This crabbing is due entirely to the 
pressure fcrces modeled in chapter II. This heading 
difference is found by computing the difference in the 
perpendicular Erejection between the FWD and AFT 
measurements and finding the arcsin of this difference 
divided ty the distance between sensors. Figure III-4 
indicates a sample of this procedure. 


The expressicn for total desired heading is given as 
follows: 

PSIDES = PSIALCC+WTSENSePSIDIF+PSIE 
where: 

PSIDIF = JY additional heading due to heading DIFference 


WISENS 


WeighTed heading difference SENSitivity gain 
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Distance Logic 
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Figure III-4 


Heading Difference Calculation 
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OO 


PSIDES = § (heading) DESired 


Thrctghcut the subroutines and main DSL programs, the 
Function CEGRAD (conversion of degrees to radians and 
tradians tc degrees) is used freely. An explanation and 


listing cf this functicn are presented in appendix A. 


"The angular velocity of the receiving ship's head is 
alsc cf ccncern in the RAS situation. This quantity may be 
thought of as similar to tachometer feedback in a simple 
serve ccntrcel system; and is necessary to damp out the 
response (the responses of this control system without this 
feedback is presented in the latter section of this 


chapter). 


The desired rudder command is a combination cf the 
desired heading, angular velocity feedback, and a rudder 
gain as fcllcws: 


Desired Rudder = (YAWD2-PSIDELC+BDOTFB) e&GN 
where: 

YAWDZ2 = heading of ship B (in degrees) 

PSILED =EFSIDES (in degrees) 

BDOTFB = VFBGeBLOT2ZD 

VFG = Velocity FeedBack Gain 


BDOTZD = angular velocity of ship B heading angle (in 
degrees/sec) 


BDOTFB = angular velocity FeedBack 


RGN = Rudder CGaiN 


The convention for rudder response dictates negative 
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rudder as kéing fright rudder, which causes positive yaw. 
This necessitates making the desired rudder the negative of 
the forcing function and feedback quantities. fThe block 
diagram cf figure III-5 presents the control loop 
measurement inputs to desired rudder command. 
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3- Optimization 


== os 


Thus far the control chcice has identified four gains 
(RSENS, &SISENS, RGN, VFBG) that must be set for f¢roper 
position attainment. The nenlinear nature of the system 
which afrears in the form of distance measurements, 
interactive forces and rudder modeling do not alicw for 
straight forward determination of these gains with normal 
Gptigal ccntrcl theory. ; 


Grossly nonlinear systems require special handling to 
determine prcper gain settings. The method chosen for this 
purpose was an optimization algorithm develored by M. J. Box 
(grcegrammed locally as subroutine POXPLX). This subroutine, 
listed in appendix A, was used to locate the cost function 
saddle pcint in four dimensional space (the dimensicns being 
the previctsly mentioned gains). The drawback asscciated 
with this method is the necessity of iterating the complete 
ncenlinear sinulaticn within function FE for every evaluation 
of the ccst function. The gains sought were found, but 
unfortunately only after 2 1/2 to 3 nours of CPU time with 
every 400 iterations allowed. 


The geéechanics involved in optimizing the chosen cost 
function include required sub-calculations in many functions 
and subroutines. Figure III-6 is. a flow chart which 


demonstrates the steps, subroutines and functions required. 


Initial optimization was accomplished for one set of 
mimeidt conditions. By locking at the RAS situaticn, a 
prekable set cf circumstances were envisioned. The scenario 


setting is the approach phase where the replenishing ship 
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Optimization Flow Chart 
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Starts. a wide approach at 0.4 ship lengths (211, feet) 
lateral displacement and 5 ship lengths (2639.0 feet) astern 
of the supply ship. The desired final position is alongside 
at 0.2.ship lengths (105.6 feet) lateral separation. Eoth 
Ships have tke same initial heading (YAW angle). The supply 
Ship is at 15 kts. (1.0 normalized speed) and the receiving 
Ship makes its approach at 22.5 kts (speed control will be 
covered later in this chapter). 


Normal ccstinc of displacement error with the integral 
of the squared error (ISE) was considered as the 
optimization tcol in subroutine BOXPLX. However, this type 
cf perfcrmance measure would weigh the initial displacement 
error equally with the final position error. This problem 
can be circumvented by comparing the displacement error toa 
pre-ccmputed reference track instead of to the desired 
displacement. For the envisioned scenario, it was conceived 
that the. ccst function should weight the distance 
displacement tEeéeavier when the ships are alongside than when 
the apprcach is started 5 ship lengths astern. 


This was accomplished by using the integral cf time 
times the aksclute error (ITAE) as the optimization 
performance measure. The reference displacement was 
considered tke desired position displacement. The object 


furcticn can then Ee written as: 


te 
CB - | te|DD-ADY| dt 


to 
wheres: 
DD = Desired Distance 
ACY = Actual Displacement in the Y directicn 


t = time 
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A perfcrmance measure that is designed to obtain good 
performance must also take into account other factors 
besides just position accuracy. Consequently, another cost 
criterion was decided upon which would also set the gains to 
reduce tke agrcunt cf rudder contrceli required when alongside. 
This particular feature is derived from the desire nct to 
over contrcl with the rudder in such clcse proximity to 
ancther vessel. The inclusion..of this. term in the 
performance measure is weighted by unity while the distance 
accuracy is weighted by a factor of 10.0. This will tend to 
allow rudder action; if the desired pesition is not 
Maintained. The final approach phase cost function for 


oktaining cptimum gains has the form: 


te 
Cra = / te(10.0¢e{DD—ADY|[+1.00«jD2j) dt 
te 
where the additicnal term is: 


D2 = rudcer response of the replenishing shir 


cc. Resuits 


In the process of deciding on the best gain definitions 
previously mentioned, many optimization runs were made. 
Each set cf gains were then simulated in a corresponding DSL 
program to cbtain performance confirmaticn. Many of these 
Tuns did not live up to expectations; causing re-evaluation 
of the ccntrcl scheme until the one presented in this thesis 
was formulated. 


Table III-1 shows the input upper and lower limits of 
Search (EU, EL), starting value guess (XS), optimum gain 
settings (Output) and associated object function value (CBJ) 
for 20.0 seccnd normalized time simulation run in function 
FE. These values were then introduced into the DSL 
Simulaticn pregram listed as program #5. The results of 
Pato simulaticn are shown in figures [II~-7 thru III-12. The 
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gecgraphic flct cf figure III-10 indicates excellent 
positioning in the lateral direction while the rudder 
resfonse cf figure III-12 shows that it settles out to a 
fairly ccnstant steady state value as the ship settles into 
its desired fesiticn. The time coordinates in all flots are 
shown in actual full scale time. | 


* 


Now that the "proper" gain settings were cbtained, more 
extensive testing of the control system was reguired. Three 
different tests were contemplated: (1), allow, a large 
perturbaticn turn of the reference ship (Supply ship), (2) 
start approach of the receiving ship (ccntrol ship) from 
G@iererent initial .conditions cf lateral and hcrizontal 
displacerzents, and (3) induce external perturbations in the 
form of wave fcrces. 


The first test was simulated by turning the reference 
Ship by ncrsgal rudder action of figure IZI-13. This -turn 
with 5 degrees rudder accounted for a total reference yaw 
change of 15 degrees. The rudder action of the ccntrolled 
ship shcwn in figure III-14 was as expected. However, the 
distance paintainment portrayed in figure III-15 was totally 
unacceptable. The maximum excursion from the desired 
distance cf 105.56 feet (0.2 normalizes distance) was §5.419 
feet Co..102 normalized distance). Variances of this 


Magnitude cannot be tolerated in the RAS environment. 


Faced with this situation, the tact chosen was to 
re-evaluate tke gains for the new scenario which is called 
the turn phase. In this phase the initial conditions assume 
Steady state positioning alongside such that the lateral 
Eesition displacement (DY) is equal to the desired distance 
Pidee56 feet (0.2 normalized) ] and that the herizontal 
Fesition displacement (DX) is O70 (alongside). 
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Scre initial ferturbation is introduced by assuminc the 
relative yaw angle when alongside is negligible. 


The same procedure was followed in obtaining gains that 
would optimize a chosen cost function. Figure III-6 still 
applies excerpt that function FEA is replaced by functicn FEB 
(listed in appendix A) to simulate the new conditions. 


Cost function criteria change in this instance since 
the ships start at the desired pesition and optimally stay 
at the same relative positions. Also, the rudder response 
to such a large turning perturbation must be free to cause 
achievement GE the desired position. Due to these 
considerations, the integral of the absolute error (IAE£) 
performance measure was chosen for the optimization 
criterion ané can te written as: 


ty 
OEJ = / JADY] dt 
t 


o 


Table III-Z shows the results of the turn phase 
optimizaticn and the comparison with the approach phase 
gains. Aqain CSL simulation was performed using the turn 
phase scenario. Figures: £fI—-16. thru | [lI-—21 portray the 
gcaphical results. The rudder response cf figure IIfI-21 
indicates very sensitive response to the interactive forces 
shewnm in figtres III-17 and Ifi-18. The lateral distance 
separaticn cf figure III-20 indicates excellent position 
Maintainment with maximum excursion error of only 2 _ feet 
(0.0038 norzalized). This minimal variaticn is well within 
that which can be tolerated in the RAS envircnment. 
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To alleviate suspicions that the response from the 
gains ottained in the approach phase could be improved by 
the gains cktained in the turn phase, a simulation of the 
approach fhase was accomplished with the new gains. Figure 
III~-22 is the graphic display of the effect of these gains 
on the approach phase lateral distance positioning. 


Gacetal analysis of the results thus far clearly 
indicate the need for an adaptive control scheme: to allow 
gain adartaticn to meet the design specifications. A _ full 
acaptive ccntrol scheme for systems of this type is outside 
the score cf this thesis. References 15 thru 22 are 
‘indicatiors cf some of the literature available for ptrsuit 


of a ccmpletely adaptive control systen. 


What was done here is develcpment of a simple algcritho 
to sense when the conditions were adequate to switch from 
cneé set cf gains te another. This may be done with the two 
sets of gains developed thus far. However, for the sake of 
Sigulaticn efficiency, a third set of gains was introduced. 
This third set amounts to a change of one approach gain 
(RSENS) which has previously been defined as the ‘range 
sensitivity cain. The simulation efficiency is increased by 
decreasing the time required for the approach phase to reach 
steady state. A consequence of this procedure is a 
Teenforcement of the need for a completely adaptive ccntrol 
scheme. 


Repeated simulation revealed that ccmmencing the turn 
(in effect switching gains), before a reasonable steady 
State was reached caused results similar te those shcwn in 
figure III-ZzZ. An increase in RSENS to a value of 4.0 when 


the lateral separation error is less than 0.05 (normalized) 
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and greater than 0.005 (normalized) forces acceptable steady 
state in approximately 1/2 the time previously reguired 
usirg a single set of approach gains. 


Subroutine S&TCH (listed in appendix A) inccrpcrates 
this simple acaptive gain schedule with a ccunter sechanisn 
to sense when steady state is reached. Further study 
indicated a need to damp the yaw oscillaticns to a greater 
extent if the yaw velocity (BDOT2D) exceeded 2.0 degrees/sec 
when the cairs are initially switched to the more sensitive 
enes of the turn phase. This is an artifical adaptive gain 
for VFBG caused by computer time restricticns prevalent in a 
full scale computer simulation where both the approach and 
turn phases are desired. If the gain switching point is 
moved up ir time, as would normally be the case in a real 


life situaticn, this damping increase would not be required. 


The results of the full scale simulation using corputer 
program #6 are shown in figures I1I-23 thru III-34. The 
aprroach phase plots. of figures III-23 thru III-2@ show 
definite improvement over that previously shown in figures 
m7 tEru Ili~—t2. Figure III-27 indicates that the 
overshcot is reduced to 10.6 feet (0.02 normalized) as 
opposed to 17.9 feet (0.034 ncrmalized) that was prevalent 
Sierigure IlI-11. 


The ttrn phase plots are shown in figures I[II-29 thru 
TiI-34 and shcw responses very similar to those shoewn 
previously in f£tgures, fil therm LLI—21. The only 
Significart differences occur in the initial responses which 
are due to the incorrect initialization when the turn phase 
Was simulated individually. 
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g.- Varying Initial Conditions 


The results oktained in the previous section are most 
gratifying tut actually incomplete. This system must work 
fer cther initial. cenditions quite different from those 
envisioned in the optimization scenario. The initial 
approach can realistically commence at points other tkan 5 
ship lengths astern and displaced by 0.4 ship lengths. 


Fy sigrulating this system with varying initial 
Eesitions, the relative efficiency and worth of the ccntrol 
system can te cbserved. This was done in successive test 
runs whose iritial conditions and corresponding plot figures 
are tabulated in table III-3. For the sake cf brevity only 
these figures required to illustrate the relative efficiency 
of the ccntrcl system are included. The correspcnding 
initial cptinizaticn simulation figures are listed for cross 
reference. The turn phase plots for all runs except 4 and 6 
exactly match that of the initial simulation and are not 
repeated kere. 


Runs 3, 5 énd 6 were accomplished to show that no 
ambiguities exist in the control scheme to prohibit adequate 
real life initial conditions. Run 3 simulates the situation 
most often encountered by this author in the RAS 
environment. This scenario starts the centrol shir dead 
astern at 5.0 ship lengths and brings it alcngside at 0.2 
Ship lengths lateral separation. 


Run 5 is a situation where the apprcaching ship is 
purposely flaced out of position on the wrong side for 
approach, the control scheme adequately corrects the 
placement error and will do so for all cases of this type, 
Prcvided that there is adequate maneuvering room astern of 


the reference vessel (in this case 2.6 ship lengths was 
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experienced which gives 1.6 ship lengths bow to. stern 
clearance). 


The purpese cof run 6 is to frovide simulation fcr an 
approach from the opposite side again disproving any ccncern 
for ambiguity in the trignometric measurement scheme 
utilized. In all runs it must be emphasized that DD is the 
positive absclute distance desired and that IS provides the 
code flag for the desired side of approach. The system will 
work with DD set to some negative guantity; but the side of 
approach will reverse itself and the positicn placement will 
be correct, Eut cn the side not desired. 


Run 4 takes the desired distance in to 0.15 ship 
lengths (€0.0 feet). This distance is usually the minimun 
desired by a prudent seaman. Again, even with this mininun 
distance, the control system performs uf to desired 
standards. TIte importance of this. run cannct be overlcoked. 
Performance cf the system at this extremum indicates that 
the gains utilized are correct for all expected conditions 
encountered in Cala seas. Figures III-59 thru III-64 


pertray tke regaining plots obtained in run 4. 


The calm sea performance of the heading control systen 
is only part of the system testing required. Of even 
greater ccncern is the adequacy of the control when sea 
State is intrcduced. Section D.2. of chapter II models the 
three ccmponents of waves with two sinusoids and a small 
Candom impulse wave. These forces were intrcduced intc the 
tctal RAS simulation as shown in computer program #7. In 
this program the wave length (WL) is set to cone ship length 
and the wave direction (WD) is -015 degrees true. This 
scenario allcews for a port turn into the prevailing séa as 


is commen fpractice in experienced RAS evolutions. BY 
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giniwmizing the perturbation forces on yaw and lateral 
direction, a smocther RAS can be accomplished thus aiding 
safety and ccmfort during the actual transfer. The wave 
force maximum is taken as 0.05685. Runs were simulated 
which used maximum wave forces in the range 0.1137 to 
0.05685, wave lengths from 0.5 to 1.5 ship lengths and wave 
directions 015 to -015 degrees off the initial replenishment 
course. The ccntrol system handled all of the perturbations 
well except fcr the cases of a wave length of 1.5. This 
length of wave with a. force of 0.05685 exceeded the ccntrol 
systems capability in that the steady state conditions were 
not met before a turn was commenced. Figure III-65 shows 
this instability in the lateral distance DY of the turn 
phase. It is felt that the modeling inadeguacies of the sea 
state develcrment of chapter II coupled with a simple 
adaptive gain scheme are the source of the problem. This 
same fphercmtenon is covered in greater detail in the 
longitudinal fesition offset testing portion of the velocity 
control section of this chapter. 


Proclems of this type also manifest themselves in some 
cases when the wave force maximum (WFMA) was close tc the 
0.1137 value. If the sea state beccmes excessive, which 
this value represents, a different gain schedule or, at 


best, a mcre ccmplex adaptive gain scheme is called for. 


The flcots produced by computer program #7 are presented 
as a representative indication of the effectiveness cf the 
ccntrol system in the presence of a sea _ state. Figure 
III~66 gives the yaw results of the approach phase which 
indicates the effect of the wave action. The correspcnding 
tudder acticn cf figure III-67 compensates to give the 
smooth lateral distance shown in figure III-68. The wave 
Beetile is shown in figure III-69 with curve 2 being WY and 
curve 3 keing WN. Curve 1 is the WX profile which was not 
used in this run but will be utilized in the speed ccntrol 
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section later in this chapter. Similar curves are portrayed 
for the turn fphase. Figure III-70 is the yaw difference 
between the two ships (remembering that the reference ship 
is nct being perturbed by the interaction forces or the wave 
forces). Figure III-71 is the lateral distance DY 
Maintained Ey the rudder response of figure III-72. The 
maximum lateral separation in the turn phase is 0.0037 ship 
Pengths (1.S5 feet). The wave profile is shown in figure 
III-73 with the same wave force curve seguence as the 


approach rhase. 


As can Fe seen from these fplots, the control system 
operates very effectively in the presence of a sea state. 
Again, the develcrpment of a much more complex adaptive gain 
scheme is reéeguired to allow exceptionally high sea state. 
It is felt that the control system presented in this thesis 
is adequate for most situaticns that are encountered in the 
RAS envircnment. Only the extreme perturbations that chance 
would allow must be accounted for in a more complex adaptive 


gain scheme. 
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B. VELOCITY CONTROL 


One advantage derived from using the linearized 
equations cf motion is the decoupling of the velccity 
compenents frcem the remaining equations cf motion. This 
allows separation of the design procedures for lateral 
separaticn ccntrol and velocity ccntrol. Section A of this 
chapter designed the lateral separation ccntrol using the 
Simple speed contrcl algorithm shown in figure III-74. fThis 
control cutgput was used directly as the ship's speed (CDCT2) 
in the mcdel simulation where no attempt was made to use the 
engine respense developed im chapter [I. Functiom SEDCTE of 
appendix A shcws the control used. 
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Figure III-74 
Non-optimum Speed Law 


Because of this decoupling assumpticn, any valid 


approach speed ccntrol can be used, if used consistantly, 
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for such a design. However, in the RAS environment, 
cciplete disasscciation is not possible. Recombination 
occurs in the interactive forces and moments which depend 
upen the longitudinal distance as well as the lateral 
distance, Consequently, speed, which is directly 
responsible fcr the longitudinal distance, has a direct 
relation to the lateral distance attainment and 
maintainsgent. 

The remaining parts of this chapter deal with the 
development cf a viakle speed ccentrol algcrithm and the 


testing cf the designed systen. 


1. Type cf Ccntrel 


== 





Whenever two ships maneuver for replenishment at sea 
(RAS), the prime ccnsiderations are the time required for 
approach and the accuracy of position keeping plus 


conservaticn of fuel. 


The necnlinear control law of figure III-75 is designed 
to maintain a preselected approach speed for mnininuo 
approach time. The proper location of the switching point 
increases the complexity of the solution since the tise of 
Switching frem this speed is determined by the dynamics of 
the neonlinear fosition attainment loop, Once this pcsition 
is reached, the speed controller is switched dawn to a 
meneart pcrticn of the control law to allow control for 
perturbations about the operating position. However, small 
perturbaticns about this operating point can be tolerated 
and, in fact, are desired to allow for conservation of fuel. 
-Selection of this dead zone is wholly dependent on the 
accuracy required fore Linal “posrtron. Frqure’ ! LLi-75 
indicates a dead zone extending to +0.001 normalized 
Bwetance which in this case translates to +£0.53 feet. 


Systems fcr which fuel considerations are net a motivating 
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factor say te designed without this part of the contrcl law 
ta allow finer tracking in the position loop. 


SPEED 





Figure: III-75 
Speed Control Law 


The speed control law as explained abcve is shcwn in 
figure III-75 for an initial approach speed cf SPDO2 and a 
final estimated reference speed of SPDO1, with ADX being the 
dynamic fosition feedback defined as the longitudinal 
distance between centers of the ships referenced to the 
controlled ship's heading. Analytically, the linear pcertion 
of the contrcl law is written as: 


SEDCTH = -ADXe(SPDC2-SEDO1) +SPDOT 


Smeectic continuation of the centrol law accounts for 
Operation cn both sides of the operating point. 
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2. Cptimization 


Using this such simplified model of chapter IIT and the 
basic ccntrcl law of figure III-75, the desired switching 
curve can be established. An optimization subroutine such 
as Subrcutine BOXPLX can be used to iteratively oktain the 
optimum switching position (S#). for representative initial 
approach speeds., Figure III-76 is a flow chart cf the 
subroutines and functions required for speed ecatroL 
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Figure. III~-76 
Optimization Flow Chart 


optimization. The major merit of this nocnlinear contrecl law 
stems frcez the predetermination of the switching pcint for 


@aa possible conditions of initial speeds. This apriori 








knowlege allows for offline computation cf the switching 
pesiticn prior to commencing the approach. The cost 
function used for optimization is the ITAE which 
accomplishes two objectives. First, it forces the approach 
tc be acccmplished in minimum time. Secondly, it insures 
that the fuel expenditure will be optimized in the 
eliminaticn cf most overshoot and bang-bang control in the 
dead zone fortion of the control law. The final value of 
the positicn error must be within the specified dead zone 
and the tergzinal speed must match the reference speed 
(SPDO1). The cost function has the following forn: 


ty 
sf (te {ADX|) dt 


te 

Table III-4 is a comparison of the optimization runs 
with various initial speeds. The values shown for Si& must 
be multiplied by the speed differential (SPDO2-SPDOt) to 
obtain the corresponding value of ADX. The max/min values 
show the -tand cf values which preduce the optimum cost. 
This range cf values is attributed to the integration step 
size used in the optimization program. Experience with this 
particular optimization program indicates that errcneous 
values of the switching point are found if the step size is 
net carefully chosen. The step size may be adequate for 
integraticn, but act for location of the switching foint. 


The fpeints obtained from the optimization runs are 
Fletted in figure III-77. These points define the ncnlinear 
Switching curve which must be stored in the computer to 
mote Oftinal oaperation of the speed ccntrol for all 
approach speeds. From here there are many frocedure oftions 
Open. These cptions have as a goal some usable forg for 
Peedicting the optimal switching point for any set of 
initial ccnditions. One may choose linear straight line 
segments with an interpolation routine, cr a closed form 


Switching curve pclynomial. Due to the availability of a 
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hybrid ccnfigured XDS 9300 digital computer and AGT-10 
qraphics terminal, the latter course was taken. 


A pclyncmial curve fitting algorithm was used to cbtain 
the reguired pclynomial coefficients of best fit. This was 
done for epelynomials of order 1 thru 5. The coefficients 
and the sum cf the squares of deviation frem the original 
peints are tabulated in table III-5. The selection cf the 
order to te used is highly defendent on the degree of 
accuracy required. In the RAS problem, the average error 
maeroduced for a first order fit is 8.0 feet(1l.07 sec), 
while the fifth order fit introduces an average error of 
1.35 feet(C.180 sec). Prior acceptance of errors intrcduced 
by an integraticn (and problem) step size cf 0.8 sec allows 
for use cf a second order fit without any degradaticn of 
Simulaticn accuracy [second order average error is 2.848 
meet (0.38 sec) J]. The graphic display of figure IIi-78 
indicates very little difference in the switching curves for 
Second tc fifth order polynomial fits. Fer the sake of 
accuracy, and owing to the computer control methods of this 
thesis, tke fifth order polynomial fit shown separately in 
figure III-79 is used for determination of the switching 


point locaticn. 





A tree test of the control law is accomplished when it 
is introduced in a computer program fcr a complete RAS 
Simulaticn. Considering the performance of this ccntroller 
in a complex envircnment of full scale RAS simulaticn allows 
maximum verification of the controller design. 


The scenario for this simulation initially fositions 
the ships such that the ship being controlled starts an 
approach 5 ship lengths (2639 feet) behind the reference 
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ship and displaced 0.4 ship lengths (211 feet) to the right. 
The desired final fposition is alongside and displaced 0.2 
ship lengths (106 feet). The heading contrcel system used is 
developed in section A of this chapter. 


The approach phase is accomplished with the speed 
desired and speed acquired shown in figure III-80 with the 
corresponding position attainment exhibited in figure 
Tet . These flots show excellent switching and oftimal 


position attainment. 


The next step is to insure that the fosition keeping 
locp will maintain the desired position with an induced 
perturbation. This is accomplished by turning the reference 
ship away frem the control ship a total of 15 degrees to 
observe the reaction of the speed control locp. The 
reference ship's turn causes the relative mction between the 
ships to cre altered, making the control ship lag the desired 
position. The nonlinear control system is designed to 
correct this situation as soon as the actual position is 
outside the limits of the dead zone. Figure III-82 displays 
the desired speed and acquired speed for the control ship. 
Figure III-83 indicates that the correspcnding fosition 
deviates frcem the desired by 0.0154 ship lengths (8.13 feet) 
at the maximum excursion. This is well within the limits of 
acceptability for such a drastic perturbaticn. 


The introduction of velocity control was accomplished 
by combining the simplified engine response of chapter II 
and the speed ccntrol law developed here. By setting the 
speed desired (SPDDES) equal to the output of Function 
SPEREC and scaling the speed error (SPLERR) £O the 
nondimensional equations of motion, the velocity lcop is 
initiated. The auxillary equations added to those presented 
in chapter II are: 
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SPCDES = SPDREC (ADX,SPDO1,SPD02,Sh) 
CDCT2 = INTGRL(U02,SPDERR¥*LUC) 


These equaticns are introduced in computer program #8 to 
peaduce figures ITiiI-80 thru III-83. 


Further system study indicates that the reference ship 
speed must be known to a fairly high degree of accuracy. 
Without apriori knewlege of the reference ship speed, a 
constant bias is introduced. The amount of bias allcwable 
defines the permissible uncertainty in the reference ship's 
initial speed. This bias can amount to as much as 0.1 ship 
lengths (€4.48 feet) for a reference speed inaccuracy cf 2.5 
kncets (0.1 normalized speed). However, it is felt that the 
reference speed in any practical situation will initially be 
Mnecwn to within 0.5 knots (0.02 normalized speed). This 
mcre practical error will introduce a bias of only 16 feet. 


Cther feedback parameters can be used to offset the 
lack of apricri kncwlege of the reference ship speed. Since 
the reference ship is tracked with a high accuracy rance and 
bearing device and the controlled ship's speed is’ measured, 
a decoupled multivariable scheme is used tc further refine 
the reference ship speed. With high resclution devices 
presently available{'*], it is estimated that this can be 
dene practically to within 0.05 knots (0.002 normalized 
speed). This would bring the offset bias to 1.6 feet; well 
within previously defined errors introduced by integration 
step size. 
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Thrcughcut the development of the heading contrcl and 
speed contrcl, the scenario has followed the condition that 
tke final position would be longitudinally alongside. 
Although this is a good assumption for ships of the same 
type, it dces not account for RAS station differences for 
different ship types. To alleviate this disparity, function 
SPDREC was redesigned to allow pre-planned cffset condition 
to exist. Function SPDOFC cf appendix A is a result of this 
redesign. 


Simulaticn runs, with a change of the speed ccntrol 
function only, resulted in some unstable ccnditions existing 
in the headirg ccentrol loop. The cause of this pheéencmenon 
stems back to the adaptive gain scheme used and the changes 
made to fcrce the control loop tc a steady state value prior 
Poe a turn. By using a favorite ploy of experienced ccnning 
officers, this problem is alleviated. The ploy is to take 
the ship alcngside and then either drop back to station or 
surge forward to station. This method is accomplished by 
setting the initial offset (XOFS) to 0.0. The final desired 
offset (XCFSD) is stored and not used until the ship is 
settled out alongside. It is subsequently used as shcwn in 
the follcwing Fortran code: 


IF (ATIME.GT.450.0) XOFS = XOFSD 


This methcd solved the gain transition problem. it did 
not, however, give a completely stable simulation run. 
Unstakle ccnditions still existed at the end of the turn 
phase. this is not suprising, considering the heading 
ccntrcol optimization method used. The set of gains 
previously fcund were for the alongside scenario only. 





Different interactive forces and moments at the cffset 
position cause these gains to be no longer cptimal-. 


By relaxing the control loop in the heading velocity 
feedback gain (VFBG), sub-optimal control at all practical 
orfset pcsiticns is achieved. The gain VFBG was changed 
frem 0.064026 to O.T in the turn phase adaptive gain 
schedule without significant less cf centrol efficiency for 
alcngside operaticn (2.3 feet maximum excursion vice 2.0 
feet previously obtained). Subroutine SWICHF of appendix A. 
reflects the gain change and offset calculations reguired. 
.Ccmputer program #9 incorporates the changes required for 
offset simulation. Table III-6 is a cross reference Listing 
of the plects obtained. From these figures, the effect of 
different Icngitudinal positicns is readily apparent. An 
offset cf 0.1, eguating to 52.8 feet, causes greater lateral 
excursions when astern (XOFSD = -0.1) of the alongside 
position than when ahead (XOFSD = 0.1). The Longitudinal 
position wmaintainment, however, is essentially the same in 
all cases. _ 


XOFSD. 


LateraL 

Distance DY 

Yaw 

Difference 

Speed ® 

Longitudinal oi ~ 

eee epee pelea 


bo = 
Note: These plots are the same as those obtained from 
computer program #8 and are not repeated here. 







Approach Phase Plots Turn Phase Plots 
a c A B c 












Table III-6 
Position Offset Testing Cross Reference 
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An alternative to the method shown here is again a 
completely acaptive gain scheme which would achieve oftimal 
control instead of the sub-optimal control settled for here. 
The alternative may become even more important if the 
ncnlinear terms of the equations of moticn are considered. 
This would couple the heading and speed control designs to a 
larger extent than encountered in the interactive forces and 


mcpents. 
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5: Wave Effects on Velocity Control 


eS 


The final testing procedure involves validation of - the 
speed ccntrcl “system in the presence of waves. This 
-perturbaticn testing continues that started in section A. of 
this chapter for heading ccntrol. In chapter II the Wx 
fcrce was acdeled thru the intermediate force IF32 as: 


* 


IFiZ = KC1eD2+NC2+KC1¢HX 


By intreducing the force in this way, a severe limitaticn is 
placed cn tke magnitude of the force. In the mariner nodel 
used, the KC1 coefficient (XDELR) is considered negligable 
cr, at best, cnly 0.00005. This translates, in the original 
equations cf motion, to a maximum speed perturbation of only 
0.C355 kts. for the wave amplitude chosen. The second 
drawkack cf this method, with even greater ccnseguences, is 
that the perturbation is introduced before the contrcl loop. 
Delay of the wave perturbation is produced making it out of 
phase with the cther wave force (WY) and mcment (WN). 








Figure III-102 
Block Diagram of Wave Introduction in Speed Loop 


mictdes ta bring about uniform introduction cf this 


Wave force, its effect is inserted just past the integrator 





of the speed contrcl loop as shown in figure III-102. This 
is coded in the DSL simulation program as: 


CDOT2 = INTGRL(U02,SPDERR*LUC) +KS2*WX 


A value of -1.0 for KS2 will give a maximum wave 
Bereturbation of Q.85275 kts. (a much more realistic 
perturbation for the high sea state simulated). Figure 
III-103 portrays the speed desired and speed acquired for 
the apprcach phase in the presence of sea state. From this 
it can be seen that the speed acquired is very dependent 
upon the sea state present. The contrcl law, however, 
presents a very stable reference for the speed locp which 
gives an approach longitudinal position (DX) plot 
indistinguishabrle from that of figure III-81. More 
prominent perturbaticn results are evident in the turn phase 
peer= of figtres LII~104 and IIi-105.. The speed response of 
figure EIZ-104 allowS a maxinum longitudinal position 
excursion of 9.5 feet (0.018 normalized position). as 
ccppared with 8.28€ feet (0.0157 normalized) in calm séa. 


These results show that the speed control system is 


very stable and corrects well for. Jarge. external: 
perturbaticns. 
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IV. CONCLUSIONS AND RECOMMENCATIONS 


A. CONCLUSICNS 


The results of this design study have been most 
gratifyirg. The Fasic concepts initially perceived for the 
RAS contrcl have been realized. The decoupled ship ccntrol 
in the RAS environment is a viable and plausible idea. This 
thesis ccntains a workable system for implementaticn of 
ccmuputer ccntrolled RAS. The achievement of 2.3 foot 
Maximum excursion for lateral distance while both ships are 
eeea. turn, and longitudinally offset by 53 feet is a 
Phenomenal achievement. Having this kind of accuracy in RAS 
cperations, can vastly increase the safety of this 


ccmplicated and dangerous maneuver. 


The approach phase of RAS can be a very hair raising 
experience. Night replenishment and sea state complicate 
the "seaman's eye" method now employed in the fleet. Having 
a system that automatically handles the apprcach regardless 
of the adverse ccnditions can, again, do nothing more than 


increase the safety of the RAS maneuver. 


Schemes for ccmputer control of nonlinear systems and 
the purpcseful introduction of nonlinear ccntrol laws are 
beccming mere practical with the technological advances in 
Micro precessors. The ever increasing number of U.S. Navy 
ships with computer systems installed, makes digital 
computer shir control realizable in the present time frame. 
A gcod wicre computer or an existing installed computer 
(Such as cne used for the NTDS system) can be used in this 
vein. Erocurement of the hardware required for this RAS 
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system can be dissipated over time periods contingent cn the 
funding available. The supply ship requires only two 
reflectors fcr the range and bearing devices stationed on 
the receiving ship. All ships can be outfitted with such 
reflectors at a minimal cost, while the bulk of the hardware 
can be introduced to the ships at regularly scheduled yard 
periods. 


In tke initial conception of this thesis, a secticn was 
planned fcr cpen ocean maneuvering. After scme research on 
this facet: of ship control, it was determined that werk in 
this area has already been documented[ 2* ][ 25 ]. The 
existence of NIDS outputs for station attainment and single 


ship contrcl systems, made design in this area redundant. 


The concept cf integrated centralized ship contrcl has 
been in the tackgrcund for over a decadef{2®]. Although 
mawen. a  icw. priority due to funding considerations, its 
implementaticn seems to be just arcund the corner 27]. 
Hamwever, a .rceview of ref. 27 indicates a lack cf RAS 
Capability. Whether this is an oversight in the article or 
neglected in the design criteria is unknown. If it has been 
neglected in the design, a very real problem has’ tceen 
overlooked. The recent incidences of ship collisions while 
conducting RAS [28] emphasizes the need for inclusicn of 
this very c¢angerous maneuver in the “Integrated Fridge 
System." Lack of technology can no longer be used as an 
excuse. This thesis and other research reports [29] have 
advanced the irplementation feasibility tc a level that 
cannot be ignored. With these projects finalized in 
practical terms, their incorporation into fleet use is the 
next imperative step. 


A major effort in this area must be made. The ever 


increasing ccmplexities of today's naval ships and the loads 
being placed cn the officers and men are such that computer 
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centrol must be used; and used now! We cannot afford the 
luxury of time to prove these system u orth, but must make 
ccncerted effcrts to get them implemented before the lives 
of 300+ men are lost. 


Whenever a complicated system such as a ship in the RAS 
Situaticr is enccuntered, many facets have tc be 
concurrently analyzed. This fact has caused inclusion of 
many diagrams in this thesis to illustrate the total 
picture. Each run, with a different ccndition, requires 
Many flots te analyze the differences in the responses and 
the causes cf the differences. The computer programs shown 
de net reflect the actual run times in the JCL shown. As 
many as twenty flots were output in these programs in the 
times listed. Analysis of the actual computation times show 
that the algcrithms run considerably under the time reguired 
for real time operation. The sampled data rate used in the 
Simulaticns was 0.11 seconds. Titsatis ywell.within, the 
realizable data rates available in even the sicwest of 
today's ccmputers and microprocessors. The thrust of this 
consideraticn is that there are ne prcblems envisioned in 
ccnverting RAS simulation to real world RAS control. 


BEB. RECCMMENCTATICNS 


In the heading control design section of this thesis, 
the need for a completely adaptive gain scheme was cited. 
Again in the velocity control section, when a longitudinal 
offset was introduced, this need became even more evident. 
The first and most important recommendation for further 
study is the development of just such an adaptive gain 
scheme. 


The linear equations of motion should te replaced with 
nonlinear equations to validate the control designs advanced 
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im this thesis. Along with this, the hydredynamic 
coefficients for the Navy's modern ships are required to be 
akle to design these control systems for todays vessels. 


It is further recommended that a ccncerted effcrt be 
made to cktain data on the interactive forces and mcments 
between shirs of dissimilar tyres and sizes. These forces 
and moments gust also be available for sea state conditions. 
In fact, tke whcle area on Sea State effects cn the various 
Seip types in the RAS situation and in open ccean 
Maneuvering needs attention. Not enough data was available 
fcr this researcher to be able to pinpoint sea state effects 
Om ship hulls. Since replenishment at sea is rarely 
ccnducted in the sterile condition of calm sea, these 
consideraticns are of utmost importance to allow testing of 


any contrcl system in the Simulation stage of development. 


The intent here is not to imply that the ccatrol 
systems pertrayed in this thesis are the best for the RAS 
scenaric, but that the procedures used can ke applied to any 
control scheme desired and bench marked to the ones 
certained here. AS previously mentioned, much meaningful 
research and éesign must be accomplished tc allow system 
reliability and, mcre important, system acceptability Ey the 
officers and men who will ultimately trust their lives to 
it. This is a task that must not be taken lightly. 
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APPENDIX A 


Due to the lengthy nature of the computer prcgrams 
presented in this thesis, many functions and subroutines 
were develcred tc simplify their presentation. This 
appendix lists these functions and * subroutines in 
alphabetical crder. The computer programs reference this 
appendix and indicate the placement of the required 
functions and subrcutines. 


A kErief description for each listing is given to aid 
the reader in determining their purpose and use. The 
following is a listing of the functions and subroutines 


contained in this appendix in the order presented: 


SCERCUTINE BOXPLX 
FUNCTION DEGRAD 

FUNCTION DELAY 

FUNCTICN FE - RUN A (FEA) 
FUNCTION FE - RUN E (FEB) 
FUNCTICN FE - RUN C (FEC) 
SUERCUTINE HDGRAS 
FUNCTION KE 

MAIN PROGRAM FOR FUNCTION MINIMIZATIONS (MINIEXPX) 
SUEROUTINE RBMEAS 
FUNCTION RKLDEQ 
SUEROUTINE SLOPES 
FUNCTION SPINIT 

FUNCTION SPDCTR 

FUNCTION SBDOFC 

FUNCTION SPDREC 

FUXCTION SWCL 

SUERCUTINE SWITCH 





SUEROUTINE SWTCHF 
SUERCUTINE TRANS 
FUNCTION XLIMIT 
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SUBROUTINE EOXPLX 


This sukrcutine was used for all optimization runs in 
heading ccntrcl and speed control. It was programmed 
lecally ard is part of the IBM 360 SSP library at the Naval 
Postgraduate Schocl. A full explaination and description is 
shcwn in the first few pages of the subroutine listing. 
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FUNCTION DEGRAD 


This function was programmed to ccnvert degrees to 
radians and radians to degrees. Ae “third: ‘purpoesie, is: tc 
ccnvert degrees to a range of 0-360. It is used extensively 
thrceughout tke programs, functions, and subroutines’ listed 
in this thesis. 
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DECUVD CMOKELT C*M2HLr et envca) 


FUNCTION DELAY 


This function is used as the time delay in the speed 
centrcl cptimization runs. It was designed to be used as an 
equivalence to DELY in DSL Sinulation.The following 
variables are defired: 


E (I) is the storage array (should be initialized 
before the first function call) 


K is the delay step count 
SPDLES is the variable to te delayed for K sters 
P is the flag for delay or no delay 


P20.0 delay SPDDES 
P<0.0 function cutput equal to SPDDES 


The functicn stores the input value (SPDDES) in E({M) and 
decrements the value in array E(I) at each call of the 
munGtIOn until the value is in the positicn of E(1). The 


value is then cutput from the function delayed K intervals. 
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FUNCTION DELAY 
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PUNCTTON FE = RUN. XK 
FEA 


This function is the simulation for heading ccntrel 
optimizaticn of the approach phase. Et, is. called by 
subroutine EOXPLX. The integration step size is 0.04 with a 
mifal tise cf 20.0. In this function all initial conditions 
are set to zero except initial geographic locaticn and 
speed. The reference ship maintains a straight course and 
the contrcl ship starts its approach 5 ship lengths astern 
and 0.4 shir lengths laterally displaced to starboard cf the 
reference shir. 


The function is referred to as functicn FEA in the 
text. 
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FUNCTION FE - RUN B 
FEB 


This function is the simulation for heading ccntrol 
Ooptimizaticn of the turn phase. It is called by subrcutine 
BOXPLX. The integration step size is 0.04 with a final time 


Ore 20.0. im this function, the following initial conditions 
are non-Zerc: 


contrcl ship rudder angle D2D & Y(14) = 8.7 degrees 

lateral displacement Y(11) = 0.2 

reference ship's speed 001 & Y(4) = 1.0 

contrcl ship*s speed 002 & CDCT2 = 1.5 (after first 
step beccmes 1.0) 


The reference ship's rudder is activated tc 5.0 degrees 


between time 4.0 and 5.0. The runs were for port side 
replenishment. 


The function is referred to as function FEB in the 
text. 
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RUN B 


FUNC TION FE 
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FUNCTICN FE - RUN C 
FEC 


This function is the simplified simulation fcr speed 
control cfptimizaticn of the switching functicn SW. ft. as 
called ty subroutine BOXPLX. The function shown is for 
approach speed of 1.1 and a replenishment speed of 1.0. The 
runs were made for various realistic combinations to cbtain 


an cptimom stitching curve. 
The rur used a step size of 0.04 and a final time of 
10.0. The twc ships were run Jlinearily with only the 


longitudinal direction and motion of any ccncern. 


Tis £uretion is referred to as. functicn FEC. in the 
text. 
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SUBROUTINE HDGRAS 


This sukrcutine was programmed to calculate the desired 
heading (FSIDES) for RAS heading control. It uses the 
outputs of sukroutine RBMEAS to calculate this heading with 
gains RSENS and WISENS. The large number of outputs in the 
subroutine call statement were made for ease of DSL printed 
output fcr tracking of simulation accuracy. 


The sukroutine also incorporates a loop to avoid 
computer precision problems in the ARSIN furction. 
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FUNCTION KE 


This function is required by all optimization runs. [It 
is the function that contains constraints for subrcutine 
BOXPLX. No ccnstraints are present, conseguently function 
KE is set equal to 0. 
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MAIN PROGRAM FOR FUNCTION MINIMIZATIONS 
MINIBXPX 


This is a generalized program which calls subrcutine 
EBOXFLX. Its main purpose is input and output of the values 
reguired in the optimization runs. This is referred to as 
MINIEXPX in the text. 
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SUBROUTINE RBMEAS 


This sukroutine measures the range and bearing of the 
ferward and after stations which is required of subrcutine 
HDGRAS. This is done with trigncmetric functions as _ shown 
in chapter II. The subroutine is specifically designed to 
circumvent any ambiguities usually associated with these 


furcticns. 


It is the basis of the decoupling of the two RAS ships 
that this thesis is based. 
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FUNCTION RKLDEQ 


This function is the Runge-Kutta-Gill forth-order 
integraticn used in all optimization runs. It is prcegrammed 
Hecaily and is part of the IBM 360 SSP Ilikrary. A full 
explaination and description is shown in the first few pages 


Gf the functicn listing. 
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SUBROUTINE SLOPES 


This subroutine contains the table look-up and 
interpolaticr scheme for the interactive forces and mcments 
presented in the RAS environment. It is long and must be 
pre-compiled for mcst of the DSL simulation programs’ shown 


in this thesis. 
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FUNCTION SPINIT 


This function was designed to aid in initialization 
problems asscciated with the DSL function DELY. The effect 
is that the functicn initializes the delay lcop until it can 
ke self-suprertive. 
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FUNCTION SPINIT 
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FUNCTION SPDCTR 


This function is the speed control used during heading 
control develcpment. It is used directly as the speed of 
the contrcel ship with information cf the speed of the two 
Ships and the lcngitudinal position ADX. It contains a 
linear furcticn at +1.0 ship lengths to a dead zone of 
+0.001 centered about the alongside position (0.0). 
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FUNC TION SPDCTR 
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FUNCTION SPDOFC 


This function is identical to SPDREC except that the 
ability te offset the alengside position (0.C) is 
mucerporated. This is the speed centrol function in its 
final develcgrent forn. 
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FUNCTION SPDOFC 
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FUNCTION SPDREC 


Pars function is similar te SPDCTR except that a 
Switching function is incorporated. This is the furction 
used for cptizization of the switching function and is used 
in the velccity loop simulated in the velocity ccntrol 
section cf ckapter III. 
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FUNCTION SPDREC 
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FUNCTION SWCL 


This functicn contains the fifth crder pfelynomial curve 
fit for the crptimal switching positicn of the speed ccntrol 
1cop. Pes range §-of valwes for SPDDIF are -0.1 to F.0 
normalized speed difference between the two ships. 
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FUNCTION SWCL 
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SUBROUTINE SWITCH 


This stkroutine contains the gains and mechanisms 
reguired for the adaptive gain schedule developed in this 
thesis. It includes the optimal gains cbhtained frcem the 


heading ccntrcl optimization runs. 
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SUBROUTINE SWTCHF 


This sukroutine is identical to SWITCH except that the 
turn phase gain VFBG is relaxed to allow for cffset 
longitudinal position placement. This is the adaptive gain 
schedule in its final forn. 
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SUBROUTINE TRANS 


This stbroutine takes the lateral and longitudinal 
geographic displacements and converts’ then to actual 
displacements referenced to the contrel ship's head. This 
is done tc gain a more realistic reference for subrcutine 
BREMEAS, sukrcutine SLOPES, function SPDCTR, function SEDCFC, 
and function SEDREC. 
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SUBROUTINE TRANS 
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FUNCTION XLIMIT 


This function was developed to allow the LIMIT function 
of DSL tc be incorforated in the optimization runs. [It is a 
saturation amplifier with a gain of 1.0, and upper limit of 
Mo, and a lover limit of XLL. 
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FUNCTION XLIMIT 
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APPENDIX B 


The final fcrm of the simulation program, with all its 
subroutines and functions, is a very complex and complicated 
maze. Tc aid in following its progression, this appendix 
ccrtains a detailed block diagram of the program with eéach 
variable listed in its computer variable name. Each page 
contains a functional part of the simulation with inputs and 
outputs shcwn cross reterenced to their origin and 


destinaticn. 


The fcllowing is a list cf the blcck diagrams ccntained 
in this appendix in the order in which they appear: 


Ship A (Beference Ship) Simulation 
Ship A Heading Simulation 
Shir A Speed Simulation 
Ship B (Ccntrol Ship) Simulation 
Ship B Heading Simulation 
Shir B Speed Simulaticn 
Subrcutine RBMEAS 
Range Measurement 
Eéearing Measurement 
Subrcutine HDGRAS 
Heading Cortrol Loop 
Auxiliary Functions 
Yaw Conversion 
Coordinate Conversion 
Feedback Loop 
Rudder Mcdeling 
Ship A Rudder 
Shif B Rudder 
Wave Generator 
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Wave Direction 
Rardem Variable 
Wave Enccunter 
Wave Components 
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APPENDIX C 


Thrcughcut this thesis subroutine SLOPES has been used 
to output the interactive forces and moments between ships 
in the RAS situaticn. This subroutine, adapted from ref. 
11, does nct ccntain a complete picture of the circumstances 
envisioned. In particular, ship's speeds cther than the 15 
kt. operating point and different ship lengths are not 
accounted for. 


As stated in-chapter II, the speed modificaticn factor 
can easily Eke applied for both ships at the same speed and 
otker than 15 kts. with the following expression: 


2 
SPDP = CDOT 


Ships replerishing with different lengths can also be 
inccrporated as shcwn in ref. 1. 


Subrcutine FAMIC listed in this appendix inccrpcrates 
these twe ideas along with a better method cf determining 
the interactive forces and moments. The curves of figures 
II-11 and II-1Z were quantized every 50 feet of DX for ail 
the DY curves shown. These peints were then used in the 
NEGS XDS-¢30C digital computer and AGT-10 graphics terminal 
Pemvoutain “a family of best fit curves. The besttit 
criteria is tased cn the sum of the error squared at each 
quantized fcint {mcdified somewhat by this researcher's 
evaluation cf test fit between pcints to eliminate spikes 
and cther ancralies). The results of this curve fit process 
is summarized in tables cC-1 and G=2, ‘which includes 
Edouglaticn cf the best fit criteria. These pclynomial 
coefficients are based on the DX distance and are ccded in 
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Best Fit Best Fit (Modified) 
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Bit Order 





¥50 £5 2.1664 8° 170.345 1.749 
Y60 a5 2.1462 14 2.1703 0.2307 
¥70 15 0.41727 12 STS 0.228 
Y80 £5 0.71421 La | 0.8758e (0-195 
¥90 £5 PUPAE os 38 14 Po293 1 0.238 
Wi00'|* 16 Lolsiz3 16 LoS7 3 0.224 
e110 15 122798 13 1.3400 0.241 
¥t20 £5 O239722 13 On 77148. 30.1835 
¥130 Ls 0.54194 18 0566737770 .170 


¥150 £5 0.26589 13 0244726, 10.139...) 





t 
¥140 eS Oa77299 | 16 1.0620 0.215 
i 
cals 
: 16 0.80572 \}0.187 








N50 0.80547 |/0.187 

N60 0.63879 |.0.167 14 0.73449 |0.179 
N70 Q.557433) "05158 14 0759042" 10.160 
N80 O243632 | 0.138 14 0.43632 /|0.138 
N90 Ge 77G85: 0.184, . 12 05 333294105290 
N100 0.48934 ee 12 0.67948 Ona 2 | 
N110 0.25702. | 0.106 14 0.59247 |0.160 | 
N120 0.29538 Oi 3 3 0.62067 |0.164 
N130 0.051807] 0.047 14 0.064323 |0.053 
N140 0.13166 0.076 10 0.30835 (0.116 
N150 oO: 11837 0.072 2 0.16549" 10'.085 
= 

Table C-2 


Interactive Curve Fit Error Analysis 


3:03 








J 





Ssuproutine FANIC as YY({1) thru YY(11) and YN(1) thru YN(11). 
An interpclaticn algorithm is used to determine the forces 
and moments at DY points between the curves of each family. 
Although all the computations are based on measurements from 
Pee, ccntrcl Ship (ship #2), the interactive forces and 
moments are also computed for the reference ship (ship #1). 


Figures C-1 and C-2 are the interactive forces and 
moments cutput to show comparison to figures II-11 and 
gi-t2. Tke speed of this run was the operating point of 15 
kts. The shirs are of equal length (527.8 feet). 


Linear interpolation of the interactive curves for 
greater than 150 feet DY distance is accomplished from this 
150 foot curve te a value of 0.0 at 200 feet, It therefore 
assumes nce fcrce and moment are present outside the 200 foot 
range. All forces and moments for DY distance of less than 
50 feet are taken as that of the 50 foct curve. These two 
endpeints are by ne means exact, but will suffice until more 
detailed cata can be gathered. Another inexact endpoint is 
preduced at the curve families limits of +550 feet. At 
these points, the forces and mcments are forced to 0.0 since 
detailed data outside of these limits was not available. A 
Side effect cf this abrupt truncation will manifest itself 
in the instantanecus commencement of the fcrces and mcments 
during the approach phase run. The endpoint variations in 
scme of the curves of figures C-1 and C-2 are due in part to 
the curve fitting routine used, but mostly to the 
differences in computer precision. (curve fits were 
calculated cr a 11 digit precision XDS 9300 while the curves 
were plotted cn single precision 7 digit IBM 360/67) 


AS previcusly mentioned, the speed medificatior for 
ether than the operating point cf 15 kts. is cnly completely 
valid fcr tke situation where both ships are at the same 


speeds. Sirce this thesis considered an approach phase 
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Figure C-1 
Curve Fitted Interactive Y Forces 
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Curve Fitted Interactive N Moments 
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where the control ship enters the interactive field at a 
speed guite different than the reference ship, some 
Bodifticaticn of .- the interactive effects shculd be 
considered. However, exact relationships are not available 


to compute tke required modification factors. 


To dispel any problems with the design of the heading 
ccntrol system, the worst case speed modification factcr was 
chosen. This factor, in effect, considers that the 
interactive fcrces and moments are derived from the ccntrol 
ship. This is accomplished in subrcutine FAMIC with the 


fcllcwing fcrtran expression: 


SPLE2 = CDOT2**2 


mo stated in chapter IL, 1 is felt that it’ is more 
accurate tc consider the interactive forces and moments to 
be modified Ly the speed of the reference ship, and can he 
coded in sukrcutine FAMIC as: 


SELE2 = CDOT1I**2 


With the scenaric followed throughout this thesis, this 
expressicr wculd equate to unity throughcut the RAS 
Situation, since the reference ship is maintained at 1.0 
ncrmalized speed (15 kts.). 


‘For the sake cf error analysis, simulation of the worst 
case modificeticn is performed. This gives rise tc forces 
and moments 2.25 times what they were in the rest of this 
thesis during a portion of the approach phase when the 
HGE@atized speed of the control ship is 1.5. Figures C-3 
and C-4 show the interactive forces and mcments for the 
aperoach phase of the simulation. The comparison plots 
which aprear in chapter III as figures III-24 and I1I1-25 
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illustrate the extent of the changes. Most notable is the 
smccther cutrut cf subroutine FAMIC. This mere 
realistically portrays the interactive effects in the RAS 
envircnment. Figure cC-5 portrays the gecgraphical flot 
which ccmfares with figure ELEE-26 without speed 
modification. Although differences exist, figure c-5 
illustrates that the interactive effects speed modification 
factor fcr tke worst case does not drastically alter the 
approach phase cutcome. The heading contrcel system cGesign 
is still valid in the face of these changes. 


For reference, figures C-6 and C-7 show the interactive 
forces and scments in the turn phase as calculated by 
sukroutine FAMIC. Figure C-8 is the turn phase lateral 
distance flet precduced. It can be seen from this 
illustration that the maximum excursion is 0.0056 normalized 
distance (2.S6 feet), well within acceptable limits. 


In summary, the designed control system will accomodate 
even the worst case modificaticn of the interactive effects. 
This insersitivity to a large range of perturbations, makes 
this control system a more viable design fcr actual ship 
installaticn. 


311 








SOOIOg AK pez4TA SAAND oseyd uUANL 


9-D sanbty 


c0'6 06° DD 'e ond 





~ “wotuu & To) 


312 














sjusUOW N peqdata esAAND sseyd UsANI 
L-) ‘eanbt ga 





05 ‘8 00 '@ 


6 006 
(ee | Meier ee eet) | Serre seen | eeeeeeoree eas 


(,014 93S) SWLL 


2 Zozwzen (x19) 


SOUS: 























S}OOTIA SATFOCASAUL peTJTpOW worT Aq soueASTq Ter197eT oseud urANY, 
8-9) eanbta 

oc 6 50°6 05°C 00°G 

Ee Se ere wee) Mee eee ee seer (eames 


(,01x 93S) awit 
(ad) adNvASIA AavIS2Q 


~ 
O 

~ 
vo 





Q-wr<zuw O> (x'o) 


314 


——— 

















_—— 


SUBROUTINE FAMTIC 


DOG SOSGOSCOSSOCOSSCOGOSCSSS ST SICOSCHSSCOGOVCSOSCeCcoe 
ANOLFONOMDAOQANUNMFPOGHODQANOATONOFM ODE QANMOTOAOMmORnDQAAIM 
dd et tet tt ION UIQ IN CICILY OO.) OO SPE SES 


HD eet FE DP ed aed Dat aman teat bad tm nt es Dt Fd fe) Pt emt ed fd Bent tes mh trem et md ft Pmt el Ga) md Peet fd ed a nt dm 
SPSSPSTPTTPSSPTAPASAASISSAPPLASTATPAALSAASTASZAAZTTES 
aIqddde ado dt ti doi a tdi dd ide 
7 iy 2 By Ol a ee ee 


% ao * 
+ 3 cy * se 
+ + al rae + 
+ + > <a + 
-~ + = ris Q * 
N + a ec + = + 
za BS lu oO * ~—-~ 4 Oa 4 
~ + N “te + O < + 
o % — OO a x< 4% 
~_ + = at * Oo ww 3¢ 
za * za uw + - tf > +t 
> + Zw + = - + 
~ + ON ke + Se st + 
N & om — 4 = % 
~ + WJ us) 3 as Ww + 
>~ + za +t =~ A te 
~ + wa wu Fe FC = a + 
Cad * =C Ste % 
~ Pig pe et m + Fa) iT] + 
> Sra any VA; =) 736, 3 
~ 3 Zaz Oe Ve) tL re e 3% 
AN i Cui +e cS w= AN % 
ae * p> aor wm wz 3% 
oO + ~~ fe See en + 
a + (=) + EVE “A. 0. 3 
(a) % wz O + A Wh m= men 
o + Ce wet ss —& HW Le * 
= 4 an Pal Si 36° Eo WH ME 
t + os el nm FH WY ft & * 
Oo + Y ee =)" *3e Ee Gin + 
fe" + Vi — * we mw Deer + 
oO + LaAaAQd = te) ees % 
- + RK RKO eee ese —-OZz * 
> + UNE seaese Qe TEN lewO + 
Q %*% WU ZUNIWMY Pie say ee AS = KO * a= 
at 3¢ Lid % <4 WrryZrO + oo Oo 
> + re LL LL Za te Sas Swit 2 m™ 
x< tt ZZ~-0CO Co t+ f& =e CLOWN * ore 
eS + dnt py <T aNe %¢ a LCL jeawt + re <a 
<a + aeZooasnstst + 2A OO Lo ++ -- 
co + SaNqoOww te << + <q O lt tw + GS => 
N + +S Ha SW (aW(aWye Aan: y > > ww tt Oo a 
=a 3 WFAA eet Ce ot OO + = 
xN + AOAWI{NeeH TTA + FFL aAOWW * ~~ © 
ot tt alan 4 GmMMw + a= mais + 9S oO 
as + TIIsTOQANWw Ww +t Wicd WYTtQa + e 1 
az + © NNAAWLL Be + OWT e~HO WwW *# oO iW 
<~ tH WN ZWNNZZO0O90Y ¢ eae ZowdZ> + UN ee 
~ e+ FS LLSNee CO ese) <>  ooaga s+ un Qa 
~ + ee (elefar eal Rea +t wiles + e 
Om * wW DudtawwZ7se *F COYRWEt N +* NAN Or 
eM + Mme Die DgeDoOOWWO + DOD AaqwH Ae: +# % # eD ui 
=—_ + REMererwoae SSD eF RPOAOYrUZ x HHN o* = 
I~ * bh OOOWGSIOOCOK- * HIdMuUWoOSst + aly <n— ——-O 
US #*+ OO Z2ZZ2Z2-2Z2UuUSTSH *€ = ara me, te bt t+ >< ANS oe 
> +t Wil Og Orie IG Me ex<e 3t oon <aNx< 
uw Ft ZF JSJIIINtee>eZZSe HF OLTOWNOZFOO * A040 *KO w 
ae ae GS) O + UNMNLDOOWAZ *# OO =e BO<k <_ 
SW +* — UNE RINGE HNN + NS xt x<I— 
ee Leal kK ;- it ag — it i EAN = 
DY + = ABNX>COSNAIN I + x it tt ir 16 
OZ *+ 2 stIOOOO>>ZZz ote + eam a xg 2 
Cw * — KKEIDOOK SKK FF + woOOo—s~Num st ~~ > 
aoz + WwW ) + Qunaa ii <i —j 
am LJ t+ tJ + eal OS aed CS) We c 
Hne: % a2 + +e KK KMS IK AO 
+ * *% 
te % + = 
hd Bia * 
* % + 
* +t *% 
OUVVUOUV OG OU OOO UU OC OOOO UU OO OUOU YO OW 


S15 

















BGOCOS HG CC Ssso OO SES oeVeGooOECooooooecce SOoeceesoeoscose 
PTOOMDOVANATOOHY DO QAaIVOAIPODMH DHOANUMN FINO DHUVGAINMTOOM—ORnOSs 
FPF POMNMMNMN WMI OO DODD 00000 E EHH ODODODDODODOAO 


tek ad es Dm es feed ed Peat met med fama Pt beet Dd Dame Gamat Pd md et ead Pt fem md med Peed Pd mes md fem amd amt md Gent ems amd fd Dame eas Pet Pat red bead nd emt tan ed Dd 
= iS pe i> Dp > ep > > p> >> a> ee > > E> Si a> > ep > i Be > > > pe > > > > > > > > > > > pe 
aAddddadadtdeoad didnt di didi dd tei dd dif odd ot to id didi data 
CLUL LL ib Lh ih ib Lh Lh LLU LL UL UL ELL iLL iLL iL LiL IL hh LLL LLL LULL LL LULL ELE 


QaANOFiUIN0OM oO 

Sete sete 

NOTIN OM OOH H Ft HH HHH HE 

HHAEH HHH THKHKHEE 

3 Ht Fe Ib He HE SE HE OC KOK OE OK OK OK OK OK 
SK ME OK OK DK OK OK OK 

Wit un wa ue at 

OQANM TNO Oo 

PION STV 0 9 tt tet at emt et 

PK HK DE MK DK OK DK BK PK HK OK OK OK OK PK OK DK DK 


NS OOD Mest 
el oN eOa 
o~ omtr Oro 
+e + e@e + O08 @ 
MoO MHNOodO MOO 
<+- <I) <t+ 
+h +h #h eR 
NS Cet NX RH 
Or OtXKM OFX XM 
+O +N +s 
wo WOO Woe 
a+ N+tO +S 
e=ttii OW + e4ui | 
MA ONW MmTrw 
em SND Wee 
Ox OTDnd ovow 
spf OND NAW 
ea em O oMoO 
Or ONRA OWN 
| e@ 1 eD 1 eo 
NO N@Ge NO e 
<1 <1orx<1 oa 
%O FO] eH wv! 
mx OMXOxXMX<XOS 
Ok OK att Ot+ 
HeN PEN KS IENX 
WO WO OWwOx+ 
O+r Nee i ores 
Ww OWwOwowag 
am OM+FMO 0+ 
OM OMweO Tri 
Aan atOorooon 
“a0 TOOT O 
Cy a) eNO ort 
On OAAMOAMND 
le 1 eDMel ey 
a9 -“AOA0-O- 
<i <1 otx<! e 
HUV HWMIOMHNS 
NX NK FH HNX | 
O+ OX OKOHD 
pet +N +0 
WO WOK HUW: + 
Or +2000 
Ow DwWwoO!] OWO 

Y Dest AD+1UHO0 + 

Ww PO STON 

QO oN OF-nHOvooe 

YF AM MOM—MMA WO 

@ 20 e0Or aenm 

Lon ornaonr 
+e + e&NiNi+ of 

~ tO COYW eOOD 
N+ Or 2000+ @ 

wm Or NSTO+ OTS 

CS ax oxeN exer 

WY 

= 

(ey 

— 

| a 

<= 

~ 

Go 

us 


SN SIN ICU) 


316 


ADDON LED Vets 
eM. « en eM «6 
AoOAsSs OW} DONO 
+ @etl +-+eci+ « e+ 
MOSM MEO MooOMN 
SFA RE Ee KR EES 
DE ee Pn Sl Del | a a 
mt >< mt SE OxXs NM 
Ox XO OFX OFM 
+OrtH# ODO +--+ eNO 
WOaA+ WCO WO! 
aAtOWw DHO N+0w 
OwW+ Pf Stu SEO 
ODuT Det mstwWo 
Onowm WMO MOmMo 
KROnoO DOM ODMOO 
msPNnre WOO Hnodao 
Ore ema etTinn 
Omnre Of AME e 
1 eMODdIi ent + aetfO 
NO ot SCI1Ge NO e+ 


<1 OrxxKx +0OTRMK +0 
TO lest te OP at UT RK 
NX OK AN XO KN XOX 
Ot st OOF Mt OF SH 
F-OKO 1 FNUXKN +S AX O 
WOt OCWWOr OUWOKO 
O+tEVNFNMN4+0 1044+ 
MmILOWMOWwOWMWOwW 
MNFMODDO+FEO—+ 44 
FOU FOES OWI CU 
TOMS HOS HONDO 
DODO Mast Om OWN Two 
eOMnO «6 eOmM~m eOODW 
QrItAOOO NS AOmnN 


| eN e+ 1 et el ol © 
AONOMHOCOAHONOC 
<1 atee<X 1 of XX] ef 


HMM KH NOMH NSM 
NX | AHN | ION I 
O+DOKHOHDXKMOHOX 
+x HORI H+ HE 
UOrA |! WOH NWO O 
O+AOwW Sf +49ON HNO 
Owo+Nnowoe | Owd + 
FO +t POT $e tue sui 
NOADDS OOMHAHODM 
NOOTPOQSHODSHN SMA 
emt ofits eONM 
ONO eO-OOOMnaO+ 
-F ENHOF OND+ OONWN 
MOF 0+00O eAOr- ee 
oO | 


ex +oux 


ANOS ANG) 


3+0.142822E+02*X1-0.533075F£+02*X2-0.103692E+02*X3+0,. 3 


ox IN ex I NX e 
PH OW H Oh +t DWH Ht Ot M+ Oa + OHHtH#O 
CON KE AIK OK OX OKLA KK HOLD OID OK SOUS KX tt 
Or Ott OFH ODOHKHO OFX OFX KO 
HrFOUv FOO Ut HAH ENS + I eI HC | 
wo woo woo wodw woo woo. 
ae nN HP OH Ee NN S$ I eR OY HT me 
SOWVOWIW Mes WU TWWOWMWW LOM WU Woot 
STN HM MOH TM HHO NOH WMS NIM OS 
SOFT TE DUD SHIM OO OI UN KN OE 
PD OPTS SAID N KOM AN ~~ ODN S&T NAIS > 
ANOS 














'@ 





DOOD OCCGDQOCOAODGO COSC DOO SCSCGe SsosoosococSoososeoo0eoao0oe 
NOPNDOM ODS ANON TO OM DOODANM POG DOAOQANOPFODM ODS SAMS OMOO 
ADDDDANOVNOOQCOCOOCOCOO0 Stasis sSsANAUNNNANNNAN MO OOM LONE 

pa ee ee ee ee ee ee ee ee ee ee ee eer) 
de Dt me amg Bad ant Ft pemed Ged Fs Dat Femed fmed aed es Fmt med Bement Det Prd Geet mad end Pema fd Dad Fendt Pad eet fred pam Damn md amet mad mt et et nd bed eed Pmt het Pad Pm) reg 
SPSS SPPPSSSTSPPSTTP PAPAS PSSPIPSTSS SP PAATPTAPATTATSTTAZSTSTTS 
aAdqdaoqadddadaddeaoidt og dod do dd i ditt dd eit dt diag 
CLE LL LLL ELL EEL UE ULE LLL ULE LUE ULL 


a= MOO OMOe eOst o ee 
Hao eNo eM™MOoO«e OFTO ONN 
00 ODD OMNRHD + o@er + 8 @ 
oe +ee Feet MOON MOS 
S0 MSD NOON —+e rH Xx $+ 
+e Kee KR HEH KEI XK CRS 
st HRA HEA KN KX N< a 
<a Xt Ka OX KH OX 
+x OFM OHKO FHtO +17 
at Fete ¢4%#O WOO! Wom 
Co WOO WOHst+ A+tOW M+O 
tO O70 tou FTIE+tO AW! 
uit NW AFD NDOWO WOW 
Nw Mest —MwO ODO frO”O 
Sm Ora OemYy trAMNUYO WOOL 
NO WMO ODN SOS KING 
NO SOO eros emoes eo 
or~ enh eC OMO 2 OMO 
Wom OOr BNO e 1 &o& t et 
oN | et | eXHt#OODNO e+ NO a 
Oe NOe NOetex!1or <x!Ig 
Io KX1O KIOwx<t OHS 1H O+ 
wo! ew HO lett KOM NxXGS 
xO NXO NXOKAOFt SH Ot 
+a OFF OF MHOC NXO +-°X 
NX FNK 4NKO1LUOX SO WO 
Oot WOt WOKOWO+N+S O90 
+Q Test NEN wWwoWwW AaAWOoO 
WO swo OWaWdOr+9O NMN+ 
NE OF + OMN+00ORPUN OOw 
Su OM NOWROSTNAN BAC 
HAF NOM NOMOHONMONGO ant 
ON TNs" TO4HOn ors Or 


Ow er™ oORO0 9Ormnn OOM 
NN ONLY SCOMMO! eO e@ ! em 
om 1 eM 1 e™@ oF#H0O=310 —OM 
aS eos AaAgmOr~x<! ef <I oe 
| x< | XS) e+ St NOM +#NOM 


INOMeHNOMe NOMA | NX | 
SM NX | AIK HOH DX OX TDM 
%HOKOK DKOHOKOAIENKM HEH +H 
ax +N XMH HPHXMHOQUIOHO LOKHN 
OFX AMWOKMMUOHKIO!A+00 O©F+00 
+HAON + AOD + HOWOUWO+ HWO ! 
wo! Two) nwesrnroertcu SNe 
V+ iu + UITO+UIOOsuUr DOIN 
OWOUOMWODasHtuTKHNOAMO OAT 
AONAITOHOMNDTNOAOANMN | rom 
ANOS OMOAIMPTINO eNNN OY el Dd Dd 
MUS 28OMN COMMS DOROO OBAMLT 
Orton Odet e+ ote + BNO 
oMNete eDU +E eS OMOD « OOF eo 
On etTOn eOON 0+0+ eOvVr+ 20 
+ eOO+ 0OO+ 2OQ0OTCO 1 FOTO 1 
TFOFHSTHO+FOLTO $s eX INK OX +E 
KEN OX HN OX HOUK OH Dei N53 Om 
HH DSN H DAO H Dette NAIK SK SOA 
NX KMONKKMONXK KO OF HO CE 
Oxt OF HH OF HOU ENT! I Pee 
teeth rate Net + WOOU WEo 
WOO WOOD WOOQWRMN + +4— + | 
De lay + 1 HO tFMNOPWUIONOWW 
ed I OT gL OS SL OD st CU) 0) ed ed OT) 
TD FLD DS OP OM Om TOO OSM MI 
DOWNS OPO SAIN Kw sO Mm HO wor 
WMT HKVUD > OLNO KNOWN Yate 
SNM ASANO BANOS BANOS RM 


QGUATICNS GF N MOMENTS 


E 


eOM «e rl ore) ei\O etm Cl eek, em 
BAND ORK) COM GAt OMAN OAIN 
+ ee! + ee | ee lee tee lee 
MOON MOO Eo Meo MES MNOo 
Ket Kee KEE KH He KI OK 
thax HEM He REI RE RR 
NM NX St aX Oe XS ee XH 
OK KN OFF KX GOHtHK OF K OF KX OHM 
+OHO .P oH Pete tH HE PH 
woo! WoO WOO WOO Woe Woz 
WrOoOw T+tO OF+0 THO M+O0 WHO 
OWw+mM Dire Mbt OWie- Mi] mus | 
WMoOousoO +FOWw OoOowW BOW BMu OGu 
wCoOMNnOe OOTY BAM NOS BOwr Wag 
NDNOODW WOT AOO NFO wWFro Wom 
Ors D CON ETO FNO NOD wen 
ero emrT eo omN 20NoO eMNco 
OMtre ONS CON O0°M OM OTM 
1 ed t eM 1 eo 1 eM | et | es 
NO ot NOe NQ 0 NOe NO oe NO e 
x<+oOr x TOTXFOLXFOrXeo <+O0 
HOt HO 1 AHO MKIH OFA OF HOF 
NXOMK WKOXXNKOXNKOXNXKD WxKD 
Or at OF BKHOK SHOOK AHOHS OtKS 
PNP a a eb ea a > ar ob 
WOrD VOX OWOKOWOHOWOH WO+ 
O+n- 1 Oral aea I NFReH In +O O+90 
OWOwW SWOWDWOWPrWOWOEWe DWE 
Otte HAOtFONTF HMO FONDS AST + 
NOWOD ASTWORAWE SEW OR Seu 
COMMN OONME-OPMODOVWVADRM HoOao 
TONM OFPWOODMAUNAMOATMHOD MOO 
eDwown TOM 2DOM EMOOD 8OVO emm 
OANN OFNOMTNORAMNOGOOARN OMM 
+ eNe + eD ot ON 0O*F 0H 0+ CO + Ov 
AOND AONORONOFOMOAAOM HOT 
<1 e+ oil el as Neth Oe Ll) el >< ae <+ 0 
HOM FTWOMNHOOMHNGMHENS HOO 
NX | AX | ANX< TAIN T RIOD $ AX + 
CHDK DHE DKOHEDKOHRDAKOHD OF D 
+ ax Ht 5 ce Pe Se i ON Pa Ge i ON PS oe PG +a4x< 
WOrKQ WO*tAWOTtAWOKAWOt WO 
DeaAD® OF00F +7AOQO0 +4900 +0 3+O 
DOF NMVOIAWSI NWSE! aAwWoO MuIe 
OOerw MN+FWOnR+WOO+UOTrT OQO+ 
OWI ODIVONNUWINODUOMeU! NOW 
DOOM CIF POPSHRAMOrFONEED MNS 
Om OMS O™NONANN SHON SNO AMON 
eran etTMO eNest erryr OM ert CO 
SD0M OFVNOCKsH00R KF ON COC 
+S OM + ODS CNH COS + 8O + OOM 
NOn eo FOAM PON “DON o+tOD OOD 
NF *O0O00+ OM! Om! OO! e Fi «@ 
OATS 1 AaAOFTO+r0TFTO1NTOI1OrO ATO 
OX 1 NK ox HN oX IN OX IN OX IN er | 
+ Om tt At OO} Oat DHMHt DHn+ O 
(>< OK CU OK Nt OK CIO XS KIRK IQS 
LGOKEHO1LOHH 1 OHH 1 OHH 1LOHH 1 OX 
+ Nee | +N +20 +230 $541 +—4 
UIMOQ OW WOO HWOO NUICO 1WOoOOe 1nuUPo 
+ 10 OFF Oe THe OH 1 m+ 
a FW US OO = OO UW LI UW WI OL UI = OU Le Ls 
VIO WOO (NOD IO OU PO MT UNG OM OF O 
SEH AIF DOM) SDN DWV we DUN er OU 
ZWOANU SPM ZOMHeNTSHZVUNU SOM 
P= HV) OWN SO) OI Re SO HLA OD HOO TOY 
ANOS HON IONS INO ORSINI) ON 


S17 











——— ~ a = - == —————— = = — —— = —— 
— en - — — ~ — - ——— — —— — =e 





SQOOSCSCOHSCSCSGOOSO GOSS CSSS OSCOSSSOSSS9OCCSS88900S6S090 
QANDIL OANO™DHOANMTMOO DKNDOANMLOAOMDAOQANYIT OOH DAOAVOTINO 
FPF PSP ST TLIO OA MWA LAWALL NO OO DODD ODOR EEE EHH ODDDODOOO 
Seas St HS Set Set St St SO et eS St Set SS Se SS St St St St eS St St St RS St SS SS St Sh Se 
ee ee ee a ed eel ele eee elt) 
StSBPSPTSSSPSTPSPSPSTSPTPTSVPSTAZPAPAPAAAASZTSASASAI TIAA STATS TSS 
CELLO LLP LPL EPL hk 


emer oom eOunr eQv a-— 
Coor B00 OFMm Get OTN 
lee tee +r coe | e |! ee 
moo Seo Moe Mo MOeo 
<rer <I! <-e+r <I <tii 
ees HA HEH HM KE 
OxXxa AX XH KX SX 
OtX OXX DEX OF OHM 
Seek i oe i a Se ty 
wood woo woo wo worn 
Or0O PHO D+O OF ++O 
Wuwer Sut Ate OW CU | 
Ost ONW MMW OO nm 
DOM MmUest ~UO WO Arm 
OMN OOF OTFm NM OOD 
ann DOr OND AIN NOW 
enwr ot OD eD@D 2D et 
OMNTN OON OMM OA OND 
| eet 1 est 1 eN toe 1 eM 
NOe NOGWe NOa NO WO e 
<+Orx!1O «KI1OrTxX<+ «KIO 
% WO +—15¢ 0 | HOlettt O #01 


NXOXNXKO AXKXOAWNXOMX<O 
StatOots OF At OF HOH 
+N Katee te NK H+ OK +O 
WOrOWO+ WOT OWOXWO* 
ats NSH OHI NEO +H 
DWOwdwO swoOwdwonu.a 
at te OLN | MmN+ MONG 1 Or | 
ANswOseu FONWMONW TOU) 
NDOMAOFST ONNYOOTONOA 
NOOOs-N HOODRHDOOMAE 
eOOv e\O° 2Ornt et oO 
OVANOMNN OsNNONDAOND 
+ DO e+ of + @Nl 0+ 00+ 02 
AONOMSOWN AONGRMONRAON 
<1 el<+ce Kl ol X+ ox<+ 


H+ NOM+ NOMS NOMHtHNOX HNO 
NX | AN + tet X< | eS | IX + 
Ot OKOHOAKOHDKOH DOH OW 
mx HRI HE HK HH HOOK +OX 
WO* NWO+# AWO + AUWIO + WO+* 
D+ QOT + AON $+ RON +N +H 
TWOEe | TWS 1! DwWoa | MWOUdUIGEe 
WO + WOO + LUE + WR | Net | 
ASAWOMNDWMOODOWOMA WN eu 
ONNSOTMNNOAINEOtNaAao 
AMO AN AOTOTFMNOWIMWY 0 
eDODOD eMNOM eOON #0 end 
GNF O/FNDONOEFONNONM 
+ e+ AMN + eta+t e+ os 
NOD eD0ON e-OR eCOONNDWM 
OMi eOM+ eOM+ 2eOD! ol+ 8 
OFO 1 O-FO+0F0+0FT00FO0 
N ex 1 N ox $QNI oX HN OX | aXe 
m0 Ht OMIM OM Ott OHO DOK OR 
tt KA XS tt SX | I | RX 
H1( Ot 1 OKRH IL OHH OF CHK Ht 
aA +890 +£0OM +O Fett FN 
Snwoorweooanwoo wo ios 
I Fer OF | re eem | ae! | 
LL) US UL CO LU LU oO UI LU Oe LN LL 
MM DHNOIDNVNADNDOO fer anor 
COND) ON TEM) OM) SON INO 
MmZUwT ZORA ZNO ZT AZANNG 
St NOM SSM tA TNT SMO HOUND 
Cy et NO (NO) ANIM BO SCM 


uw tt pet et oOo 
<< me ee 
To 1b KK Oo 
— hs feet at LOL 
oO St te ~™N™ 
-O <x x -—— 
(a>) —— — — ey 
Se Ht Kt >> 
ot ted liek 
~~“ @ —~— |! CcCSsSe 
ah >~>z ee — se 
a ine oat ~2z COWW 
© om OC -_—- +> One 
- — Of—f-- Cake 14 —_— 
~ ow +--+ - —+ + 
> cond te ooo 
(=) Cc eno —_ eee 
=< °© cag ea Ones 
>~o own >> >> o-~ 
Qasr Omit —_—-_— = Nite 3 
aca | Col oe Leelee +++ fet tome 
~tHt eox xc — = ee = he mel et 
Nam hh © Cmte eet OO 
axa eer Z2t KFS KZ rset Se rs ee 
te OD eh Oe ee re > OO 
“out md Fed Pad Ged fd Pd (=| 
“ux «KKK KH HONWCHNONnHOamNNOWn 
a ———~—— | aed - -~— - 
> II Len I ce | ce coe] oe WD mee ION 


OO Lie Se SOR SZ Ke S2L re ZO 
TL OK at FR at at at st De OO SO ee Om 


N (a8) Tr uY\ Se) 


38 








SOCOOCSCSCOOOCGeooeoeo 
SVHOSANMTUOM OWVOQaANMO TWO 
DDMAAWAYMDAAARACGHOSOOOOo 
Snes SSS SH ACN IAI OIAY 
ES ed et mee fee Pe re ems et beg ed ed fet Pe mt et 
SSSSSSSTSzTSTSzszSaATTzZzTS 
aIdqadddadodcoet ei tc 
CLEL UL LLL IL ULL LLL kL LLL LLL 


LAW 
oo LAL 
td Oo 
WLU Tea 
oo UU 
ee oo 
at e eO 
+ + at ot 
NIN H#O 
fakes Hate 
oj ey ak 
<x< — ie) 
C+ + << OD 
RAN +H 
(ou folks Se 
oad aag 
oOaa QAO e 
WN AAR 
—t + NW e 
Nee HEE SNA 


ag OosZz N may ee ZZ 
ee U>> UW Heze >>>> 


Oo Wl lil ax<a> >> Het Ve Oi 
ze +! oz 

nbaesznunuZz Onvrdanwnuninae 
2_— ihe ~D = 

SNe NN I ae INS 


22WUU>2 JOY ZUw>Ksezw2! 
eH HZ KO} HH 


~ wm Qo 


39 











COMPUTER PROGRAM #1 


This program incorporates the ship dynamics of two 
identical Mariner hulls. These hulls are superimposed in 
space to zgllcw fcr comparison of the effects contributed to 
rudder mcdeling differences. [In this particular run a step 


and ramp rudder were compared in chapter II. 


Ancther kenifit of this pregram is tc set up the two 
identical shirs required for the RAS simulations in chapter 
III. Basicaliy, verification cf the models in three degrees 


cf freedem is acccmplished for the Mariner hull chosen. 


The rlcts froduced in this run are shown in figures 
mie? and TI- 3. 
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COMPUTER PROGRAM #2 


This prcegram models a practical rudder response fecra 
Mariner shir type. The rudder limits (stops) are set at +30 
degrees anc the cate of response is limited toa +2 
degrees/seccnd. A scale factor (LUC) is introduced to 
modify the respcense to match real time of the mariner hull 


chosen. 


Twelve passes thru the program are accomplished to 
ccnform to different sets of initial conditions and final 
desired rudcer conditions. The flots produced in this run 


are shown in figures II-5 and II-6. 
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COMPUTER PROGRAM #3 


This pregram models a reduced order (first order) gas 
turbine froprlision plant for an input-output relationship. 
The program does not scale tke plant to the mariner hull 
used. This was done when introduced into the main 


Simulation prcgram first listed as computer frogram #8. 


The time delay (P) is assisted in initialization by a 
dual feed into the system; one thru the delay itself and one 
directly intc SPDIN. The program can be modified to ccmpare 
a family cf curves. by introducing the following sequence 
into the TERMINAL region: 

INTGCER NUMB 

IF (NCUR.EQ.NUMB) CALL ENDEW(NPLOT) 

IF (NCUR.NE.NUME) CALL RERUN 

NCUR = NCUR + 1 
where NUME is the number of curves desired (less than or 
equal to 10) which is set with a PARAM statement. The 
comparison is dore on the conditions set in the terminal 
region [i.e. decrement or increment the system gains (€g. G 
=——G, + 0.C2) j. 


fie pect procuced, by this (run. is shcwn- as part of 
figure II~-10. 
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CCMPUTER PROGRAM #4 


This pregram models a simplified wave Simulation 
ccmpesed cf two superimposed sinuscids (fundamental and 
seccnd harmcnic) and a small random wave. The model is 
inherently scaled - to the mariner nondimensional 
characteristics. Introduction of these waves is accomplished 
in computer program #7, Subroutine DEGRAD is shcwn in 


appendix A. 
The sea state force plots in the dimensions cf the 


three degrees cf freedcem produced by this run is shown in 
figures; II—17 thru II-22. 
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COMPUTER PROGRAN #5 


This fiecgram uses the mariner hull model FIese 
introduced in computer program #1 and the control system 
designed in chapter iII to simulate the approach phase of 
BAS. The subroutines and functions that are to be inserted 
frcem appendix A can also be dore in object cede by changing 
the word FORTRAN to OBJECT and placing pre-compiled decks in 
the same locations. In fact, due to the Ilcng length of 
subroutine SI1OPES, this must be done to be able to run the 
Sipulation with the DSL default job control language (JCL) 
presently installed at the Naval Postgraduate Schocl IBM 
36C/67. 


The flets froduced by this run are shown in figures 
mes thru EIt-13. By changing the gains and introducing 
the follcwirg code, the plots of figures III-14 thru III-19 
are prodtced: 

DILES = 5.0*STEP (8.0)-5.0*STEP (9.0) 
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COMPUTER PROGRAM #6 


This pregram combines the approach and turn phases of 
ccgputer prcegram #5. The added subroutine is a result of 


Simulaticn requirements to switch between adaptive gains. 


This .ctur produced the plots of figures III-22 thru 
777-34. By substituting the initial conditions of table 
IfI-3, tkis program produced the plcts cf figures III-35 
thru TITr-€4. 
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COMPUTER EFROGRAM #7 


This pregram combines the calm sea simulaticn of 
cciputer program #€ with the wave simulaticn of computer 
program #4 tc simulate the model and control system in a sea 
state. The waves are introduced thru the rudder 
ncndimensicnalized ccefficients as shown in chapter II. 


The flcts produced are shown in figures IITI-66 thru 


III-73. Figure II1-65 was produced with the same program by 
setting Wl=1.5. 
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COMPUTER PROGRAM #8 


This pregram incorporates a fifth ocrder pcelynonial 
curve fit speed ccntrol switching functicn te give oftimal 
lcngitudinal fositioning. The scenario is the same that was 
used in the Gesign of the heading control development. The 
low order sucdel of the gas turbine propulsion plant was 
used. 


The flcts produced are shown in figures III-8C thru 
Prt —o 3. 
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COMPUTER PROGRAM #9 


This pregram introduces a longitudinal pesiticn cffset 
capability. The sethod takes the control ship tc the 
alcngside ecsition until 450 seconds into the run. After 
that time, with the ship steadied, the cffset position 
desired (XCFSD) is switched to the desired offset. This 
method negates some of the transient oscillations which 
cause unstacle conditions in the apprcach phase. A 
seccndary change is the use of subroutine SWTCHF instead of 
Sukroutine ShTCH developed in the heading control section. 
This new subrcutine relaxes the heading velocity feedback 
gain (VFEG) tc allecw turn stability in the turn phase. 


The plots produced by this prcegram are shown in figures 
pero thru TIT10'1. 
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COMPUTER FROGRAM #10 


This pregram incorporates the sea state first 
prcgrammed in computer programs #4 and #7. The WX wave 
fcrce, hcwever, is introduced at the end of the propulsion 
loop. to allcw more realistic perturbations. This is the 
final fcrm cf the complete heading and speed ccntrol 
systems. Tc cun this without a sea state, set WFMA tc 0.0. 


This frcegram produced the plots shewn in figures 
mer—103 thru ITrI-105. 
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